FLUID  DYNAMICS  OF  CIRCUMSTELLAR  MATERIAL 
ASSOCIATED  WITH  Be  STARS 


By 
Thomas  Harlow  Morgan 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE  COUNCIL 

OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT  OF 

THE  REQUIREMENTS  FOR  THE  DEGREE  OF 

DOCTOR  OF  PHILOSOPHY 

UNIVERSITY  OF  FLORIDA 
1972 


To  KYC  and  HEN — both  of  whom  pushed. 


ACKNOWLEDGEMENTS 

I  am  debtor  both  to  the  Greeks,  and  to  the  barbarians; 

both  to  the  wise,  and  the  unwise. 

Romans,  il4 

I  would  like  to  thank  my  parents  for  their 
encouragement  throughout  my  education. 

I  am  indebted  to  Dr.  D.  C.  Swanson  for  first  interesting 
me  in  physics  as  well  as  to  Drs.  S.  S.  Ballard,  C.  F.  Hooper, 
Jr.,  and  F.  E.  Dunnam  for  help  and  advice  in  the  summer  of 
1970.   Dr.  F.  B.  Wood  introduced  me  to  variable  stars,  and 
Dr.  K-Y  Chen  suggested  the  topic  of  this  dissertation. 

Mr.  W.  W.  Richardson  drew  the  figures  that  appear 
here.  The  members  of  my  committee,  particularly  Drs.  K-Y 
Chen  and  A.  G.  Smith,  read  and  criticised  the  rough  draft. 
Dr.  R.  S.  Leacock  also  read  sections  of  the  rough  draft, 
and  in  addition,  offered  many  practical  suggestions.  Mr. 
H.  E.  Nutter  was  of  great  help.   I  am  indebted  to  them  all. 


•  •  • 

iii 


TABLE  OF  CONTENTS 

Page 

ACKNOWLEDGEMENTS *■**• 

LIST  OF  TABLES ix 

LIST  OF  FIGURES X 

KEY  TO  SYMBOLS xii 

ABSTRACT XV 

CHAPTER 

I      INTRODUCTION  1 

Discovery  and  Taxonomy  of  Emission-Line 

B  Stars.  . 1 

A  Brief  History 1 

Nomenclature  2 

Characteristics  of  Be  Stars 3 

Spectral  and  Luminosity  Class ....  3 
Galactic  Distribution,  Membership 
in  Multiple  Star  Systems  and 

Incidence  Among  Early-Type  Stars  .  3 

Rotational  Velocity 4 

Spectral  Variations 5 

The  Struve  Model  6 

Origin 6 

Line  Profiles 10 

Hydrodynamical  Models  for  the  Flow  .  10 
Spectral  Variations  and  their 

Possible  Causes 11 

Modifications  of  the  Struve  Model.  .  12 

Major  Topics  of  this  Work 13 

II    HYDRODYNAMICAL  APPROACHES  AND  STEADY-STATE 

SOLUTIONS 15 

Background 15 

iv 


CHAPTER  Pa9e 

The  Hydrodynamical  Equations  and  Notations.  .  .  15 

Static  Solutions .  18 

The  More  General  Solutions 20 

Form  of  the  Circular  Velocity  Law 23 

Final  Comments 27 

III     APPLICATION  OF  HYDRODYNAMICS  TO  THE  ENVELOPES 

OF  Be  STARS 29 

Disk  Dimensions 29 

Theoretical  Foundations  of  Hydrodynamics.  ...  31 

Validity  of  the  Boltzmann  Transport 

Equation 31 

Relaxation  Times 32 

The  Hydrodynamical  Prescription  35 

The  Hydrodynamical  Equations 35 

Navier-Stokes  Equations  and  Similarity 

Numbers 37 

Importance  of  Viscous  Term  in  the  Equation 

Of  Motion 41 


IV     MATHEMATICAL  ANALYSIS  OF  THE  HYDRODYNAMICAL 

EQUATIONS 44 

Preliminary  Remarks  44 

Linearization  of  the  Equations  for  Temporal 
and  Angular  Dependence 45 

Introduction  of  Time  and  Angle 

Dependent  Terms 45 

Selection  of  a  Subscript- Zero 

Solution 50 

Considerations  Concerning  Approximations 
Near  the  Plane 53 

Additional  Relations 56 

Integration  Over  z 58 

Further  Restrictions  on  the  Flow.  ...   58 
Relations  Among  the  Flow  Variables.  .  .   64 

v 


CHAPTER  Pa?e 

A  Reprise. 68 

V     COMPUTATIONS  AND  RESULTS 70 

Introduction  to  the  Computations  ....  70 

The  Equation  for  ic 71 

The  Consequences  of  Imaginary  Roots.  71 

Initial  Values  for  the  Parameters.  .  72 

Machine  Computation  of  the  Roots  .  .  74 

The  Case  n  =  0 75 

Distinctive  Aspects  of  the  n  =  0 

Case 75 

Variation  with  w,r  and  A 75 

The  Case  n  =  1 78 

Preliminary  Remarks 78 

The  Effect  of  c, 79 

The  Effect  of  AT 93 

Extreme  Cases 93 

The  Case  n>l    107 

General  Features  107 

Large  n  and  the  Effects  of  c1  and  r.  107 

Sample  Flow  Variable  Solutions  114 

Comparison  with  Observations  115 

Actual  Stellar  Values  for  Vs  .  .  .  .  115 
Qualitative  Aspects  of  the  Predicted 

Spectral  Variations 116 

Comparison  with  Observations  Reported 

in  the  Literature 118 

VI     EQUATION  FOR  ENERGY  TRANSPORT 121 

Form  of  the  Equation  of  Energy  Transport  121 

The  Distinctive  Nature  of  the  Energy 

Equation •  121 

An  Approximate  Form 123 

Relative  Importance  of  the  Terms  of 

the  Equations I25 


VI 


CHAPTER  Page 

A  Linearized  Equation 129 

Separation  of  the  Time  and  Angle 

Dependence 129 

An  Equation  for  the  Subscript-One 

Terms 131 

Relationship  to  Previous  Work 134 

Reduction  to  Comparable  Form  ....  134 

The  Question  of  Consistency 137 

A  Reconsideration  of  the  Approximate 

Energy  Equation   138 

Non-Static  Solutions  and  the 

Adequacy  of  the  Equations 138 

Final  Comments 140 

VII   FINAL  COMMENTS 142 

Stability  of  the  Steady-State  Flows.  .  .  142 

Conclusions 144 

APPENDIX  A 148 

Exchange  of  a  Differential  and  an 

Integral  Operator 148 

APPENDIX  B 151 

Separation  of  the  Equations  and 

the  k  Relation 151 

APPENDIX  C 155 

Evaluation  of  the  Determinant  formed 
from  the  Equations  Relating  U  .  , 

Cnk'  «**,*•.".' -  -  •  •  -  155 

APPENDIX  D 157 

The  Program  CALSOL 157 

APPENDIX  E 165 

Pertinent  Integrals 165 


vi  x 


Page 

APPENDIX  F  ......  . 168 

An  Equation  for  r  __ 168 

BIBLIOGRAPHY 171 

BIOGRAPHICAL  SKETCH 174 


Vlll 


LIST  OF  TABLES 

Figure  Page 

1  Characteristic  Central  Star  and  Disk 

Values 30 

2  The  values  of  the  Coefficients  Appearing  in 

the  Energy  Equation. 128 

3  The  Program  Names  for  Important  Quantities  .  .  159 

4  Main  Program 160 

5  Subroutine  ALPl 163 

6  Subroutine  BET  1 164 


ix 


|PER0  |    for  small  values  of  c1.x 81- 

fc-Jand     . 
|PERjfor  larger  values*  of  cx 83 


LIST  OF  FIGURES 

Page 

Fxgure 

1  A  Be  star  viewed  equator-on  and  pole-on  ...    8 

17 

2  The  coordinate  system  

3  The  radial  dependence  of  JPER|  for  n  =  0        ^ 
and  for  large  c1 

The  radial  dependence  of  JPER-jJ  and 

2 
The  radial  dependence  of  | PER, | and 

'a'  ■       " 

6  The  radial  dependence  of  |PER3|for  four 
values  of  c1 

7  The  radial  dependence  IPER.J  and  |PER2| 
for  small  values  of  r  . Bb 

8  The  radial  dependence  of  |PER.J  and 
|PER2|for  large  values  of  r yu 

9  The  radial  dependence  of  |PER3|  for 
small   values  of  r  

10  The  radial  dependence  of  |PER3(  for  three 
values  of  T  

11  The  radial  dependence  of  all  three  periods 
for  large  c1   and  r   ....  

12  The  radial  dependence  of  IperJ  and|PER2| 
for  small  values  of  c1   and  r a:* 

13  The  radial  dependence  of  the  three  periods 
for  three  different  values  of  n iy)* 

14  The  quantity  IPER-J^  as  a  function  of  n  .  .  104 

15  The  effect  of  large  parameter  values  on 
the  periods  when  n  is  large  


Figure  Page 

16  The  radial  dependence  of  u  -(0°.)  and 

11^(0")  for  various  times'..  . HI 

17  The  radial  dependence  of  U., (90°)  - 

uil(270°)  for  various  tirae§7  .......     113 


XI 


KEY  TO  SYMBOLS 

In  order  to  avoid  an  excessive  number  of  symbols, 
extensive  use  was  made  of  subscripting.   The  subscripts  fall 
mainly  into  four  classes.  First,  subscripts  are  used  to 
denote  a  component  of  a  vector.   For  the  cylindrical 
coordinate  system  in  use  here  (see  Figure  2) ,  these  subscripts 
are  a),<j>,  and  z.  Second,  the  subscript  c  indicates  a 
dimensioned  quantity.  Third,  a  capital  letter  subscripted  o 
is  a  characteristic  dimensioning  quantity.  Finally,  a  subscript 
o  appearing  with  a  lower  case  letter  indicates  a  steady-state 
cylindrically  symmetric  term,  and  the  subscript  1  indicates 
a  term  with  temporal  and  angular  dependence. 

The  abbreviations  for  units  follow  the  usage  of 
the  Manual  of  Style  for  the  Astrophysical  Journal.   Symbols 
whose  usage  is  confined  to  a  page  or  a  section  will  not  be 
listed  here. 

Important  Symbols 
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V     velocity 
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Temporal  and  angular  variations  in  the  motion  and 
distribution  of  circumstellar  material  about  typical  Be  stars 
are  studied  by  means  of  a  simple  hydrodynamical  approach.  The 
validity  of  any  hydrodynamical  study  of  the  dynamics  of  the 
envelope  is  first  examined,  and  the  problem  is  found  to  lie 
within  the  realm  of  hydrodynamics  for  typical  disk  conditions. 
The  problem  is,  however,  close  to  the  limit  of  validity  of 
that  discipline.  Further,  Euler's  equation  may  be  used  as 
the  equation  of  motion  throughout  most  of  the  gaseous  disk 
surrounding  the  central  star. 

Each  flow  variable  is  written  as  the  sum  of  a  known 
steady-state  axially-symmetric  term  and  an  unknown  term 
containing  the  temporal  and  angular  dependence.  The  squares 
of  these  latter  quantities  are  assumed  small.  Using  this, 
the  equations  of  motion  and  continuity  are  reduced  to  a  set 
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of  linear  equations  for  the  terms  with  temporal  and  angular 
dependence.  The  known  steady-state  axially-syrametric  terms 
are  taken  from  the  literature  (Limber,  D.  N.  1964/  Ap.  £«* 
140,  1391).   This  solution  is  a  static  isothermal  one.   The 
geometry  of  the  problem  suggests  several  approximations  which 
may  be  used  in  the  linearized  equations;  these  approximations 
introduce  two  parameters.  Solutions  can  be  found  to  this 
somewhat  approximated  version  of  the  linearized  equations. 
Since  the  boundary  conditions  have  not  been  treated,  a  fully 
determined  set  of  solutions  is  not  possible.   The  dependence 
on  both  time  and  angle  is,  however,  explicit. 

The  angular  dependence  enters  the  solutions  in  an 
extremely  simple  form  involving  an  integer  n,  but  understanding 
the  temporal  dependence  requires  computation  of  the  periods 
characteristic  of  the  temporal  variation  at  a  given  location 
in  the  gaseous  envelope  and  for  a  given  choice  of  the  two 
parameters  discussed  above.  The  IBM  360/65  at  the  University 
of  Florida  was  used  to  perform  these  computations.  Under 

most  conditions  the  characteristic  periods  are  found  to  be 

3/2 
real.  Further,  their  dependence  on  n  and  uis-io   /n  where 

co  is  the  distance  from  the  center  of  the  star  in  the  equatorial 
plane.  These  calculations  predict  temporal  variation  down  to 
fractions  of  an  hour,  and  in  qualitative  agreement  with 
relevant  observations . 

The  third  equation  of  hydrodynamics,  the  energy 
equation, is  examined  insofar  as  it  can  be  in  the  absence  of 
an  effective  treatment  of  the  radiation- fluid  interaction. 
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This  approach  leads  to  inconsistencies  with  the  results 
of  the  study  of  the  equations  of  motion  and  continuity. 
It  is  concluded  that  in  the  absence  of  some  accommodation 
of  radiation  effects,  no  equation  for  energy  transport  can 
improve  the  existing  solutions. 

Finally,  the  stability  of  steady-state  flows  is 
examined.   The  steady-state  flows  are  not  absolutely 
unstable,  but  the  most  general  solutions  appear  to  be 
periodic  in  character. 
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CHAPTER  I 
INTRODUCTION 
.Discovery  and  Taxonomy  of  Emission-Line  B  Stars 

A  Brief  History 

In  August  1866  Father  Secchi  observed  HB . in  emission 
in  y  Cas.   Over  the  next  half -century  Campbell,  Frost,  Miss 
Cannon,  and  others  observed  hydrogen  in  emission  in  other 
stars  of  what  would  now  be  called  spectral  type  B  as  well 
as  in  members  of  the  remaining  early  spectral  types.  Miss 
Cannon  appears  to  have  been  the  first  to  observe  variations 
in  the  emission  features  in  one  of  these  stars  (Cannon,  1898) . 
Pickering  (1912)  published  the  first  catalogue  of  early-type 
emission-line  stars;  at  that  time  there  were  94. 

Further  work  was  in  two  directions:  more  emission- 
line  objects  were  discovered  and  individual  stars  were 
studied  in  detail.   The  catalogue  of  Merrill  and  Burwell 
(1933,  1943,  1949)  contained  over  200  B-type  emission-line 
stars.   Detailed  spectral  studies  led  Struve  (1931,  1942) 
to  interpret  the  spectral  features  characteristic  of  Be 
stars  in  terms  of  emission  and  absorption  in  a  gaseous 
disk-shaped  envelope  surrounding  a  rapidly  rotating  central 
star;  this  interpretation  has  not  changed  substantially 
since . 

In  recent  years  both  long  term  variations  (Limber,  1969) 


and  shorter  term  variations  (Lacoarret,  1965;  Slettebak, 
1969;  Hutchings,  .1969  and  1970;  Hutchings  et  al . ,   1971)  in 
the  characteristic  emission-absorption  features  have  received 
attention.   The  most  recent  listing  of  Be  stars  (Jaschek 
et  al . ,  1971)  contains  1900  objects. 


Nomenclature 

The  basis  of  modern  systems  of  spectral  classification 
is  the  system  adopted  for  the  Henry  Draper  catalogue  (Cannon 
and  Pickering,  1918-24) .   The  original  Harvard  arrangement 
of  the  spectra  was  in  order  of  decreasing  ratios  of  intensities 
of  the  hydrogen  Balmer  lines  to  the  intensities  of  a  number 
of  other  absorption  lines.   The  Harvard  observers  soon 
realized  that  a  rearrangement  of  the  original  alphabetical 
order  into  the  present  spectral"  sequence — 0,  B,  A,  F,  G, 
K,  M — represented  an  ordering  in  temperature  with  decreasing 
temperatures  from  O  to  M. 

Among  the  HD  objects  typed  B,  those  possessing 
hydrogen- emission  lines  were  marked  with  an  additional 
lower  case  p —  for  peculiar — following  the  B.  Modern 
classification  uses  the  notation  Be  for  such  stars,  that 
is,  stars  with  B-type  spectral  features  plus  emission  in 
one  or  more  of  the  Balmer  lines.   The  suffix  p  is  now 
used  only  to  indicate  the  presence  of  annomalously  strong 
or  weak  lines.   If  a  Be  star  shows  sharp  metallic  absorption 
lines  superposed  on  the  Be  spectrum,  the  star  is  called  a 
shell  star. 


Characteristics  of  Be  Stars 


Spectral  and  Luminosity  Class 

In  this  work  the  term  "Be  star"  applied  to  hydrogen 
emission-line  stars  whose  underlying  spectral  types  are 
early  of  mid  B  (B0-B7)  and  whose  luminosity  classifications 
are  either  II,  IV,  or  V.   This  is  somewhat  at  variance  with 
much  of  the  literature,  which  considers  Be  stars  to  be 
basically  main-sequence  (luminosity  class  V)  objects.   But 
Mendoza's  (1958)  comprehensive  luminosity  class  survey  of 
Be  stars  in  the  Merrill  and  Burwell  catalogue  simply  placed 
the  stars  in  two  groups — very  luminous  objects  (I,  II)  and 
"less"  luminous  objects  (III,  IV,  V) — on  the  basis  of  a 
pronounced  dichotomy  in  his  observational  material.   In 
fact,  strong  doubts  (Underhill,  1966)  have  been  expressed 
as  to  the  existence  of  a  systematic  distinction  between  MK 
luminosity  classes  IV  and  V  in  this  early  spectral  region. 
Absolute  magnitudes  of  Be  stars  determined  on  the  basis  of 
cluster  or  association  membership  appear  to  fall  about  one 
magnitude  above  the  zero  age  main  sequence.   Finally  the 
appellation  shell  will  be  similarly  restricted  although 
objects  with  A  and  even  early  F  spectral  classes  are  some- 
times called  shell  objects. 


Galactic  Distribution,  Membership  in  Multiple  Star  Systems, 
and  Incidence  Among  Early-Type  Stars 

Mendoza  (1958)  took  for  his  sample  all  the  stars  in 
the  Merrill  and  Burwell  catalogue  whose  spectral  types  fell 


in  the  range  from  BOe  to  B2..5e  and  whose  luminosity  classes 
were  III,  IV,  or  V.  He  found  that  the  Be  stars  populate  in 
a  roughly  uniform  manner  the  spiral  arms  as  delineated  by 
the  0-B  associations.   Be  stars  are  found  in  0-B  associations 
but  are  not  strongly  concentrated  in  them.   They  are  to  be 
found  both  in  galactic  clusters  (x  Persei,  Coma)  and  also 
in  binary  systems  (a  Tau,  <j)  Per,  V  367  Cyg)  . 


Rotational  Velocity 

Be  stars  are  the  most  rapidly  rotating  class  of 
stars.   The  most  recent  catalogue  of  observed  rotational 
velocities  (Bernacca  and  Perimotto,  1971)  lists  189  objects 
in  this  spectral  class  whose  rotational  velocities  have 
been  measured.  The  mean  rotational  velocity  of  the  group 
is  262  km/s,  where  the  observed  rotational  velocity  is 
the  true  rotational  velocity  times  sin  i,  i  being  the  angle 
between  the  axis  of  rotation  and  the  observer's  line  of 
sight.  On  the  assumption  that  the  axes  of  stellar  rotation 
are  directed  randomly  in  space  one  finds  the  true  mean 
rotational  velocity  (Huang  and  Struve,  1960)  to  be  344  km/s. 
Traditional  methods  of  line  profile  analysis  used  to 
determine  the  observed  rotational  velocities  are  now  believed 
(Collins,  1970;  Hardrop  and  Strittmatter,  1968)  to  underestimate 
the  actual  values  by  as  much  as  40  percent.   If  this  belief 
is  true  the  Be  stars  are  rotating  at  large  fractions  of 
their  equatorial  break-up  velocities.  Balmer   emissions 


occur  in  approximately  15  percent  of  all  B-type  spectra 
(Underbill f  1966) . 


Spectral  Variations 

Spectral  variations  occur  in  most  Be  and  shell  stars. 
Following  the  terminology  of  McLaughlin  (1961)  the  spectral 
variations  fall  into  three  major  classes: 

(1)  appearance  or  disappearance  of  a  shell 
absorption  spectrum, 

(2)  changes  in  the  ratio  of  the  intensity  of  an 
emission  line  to  the  neighboring  continuum  (E/C) , 

(3)  changes  in  the  relative  intensities  of  the  red 
and  violet  components  of  the  Balmer  emission 
lines  (V/R) . 

The  classic  example  of  the  first  variation  is 
Pleione  (28  Tau) ,  which  has  been  observed  spectrographically 
since  1888.   Pleione  was  first  classified  as  an  emission- 
line  star,  but  by  1905  it  showed  no  emission  features  and 
possessed  a  normal  B  spectrum  (with  large  rotational  velocity) 
In  1938,  emission  lines  were  present;  a  short  time  later 
the  star  displayed  a  shell  spectrum.   By  1954  the  emission 
features  were  lost.   Recent  observations  (Sharov  and  Lyutiv, 
1972)  indicate  that  Pleione  is  entering  a  new  cycle  of 
activity . 

E/C  and  V/R  variations,  though  strictly  periodic 
only  when  they  reflect  orbital  motions  of  stars  in  binary 


systems,  show  much  quasi -periodic  behavior.   The  time  scales 
for  these  variations  may  be  on  the  order  of  a  few  years  or 
as  short  as  days  or  fractions  of  a  day  (Slettebak,  1969; 
Hutchings,  1969  and  1970;  Peters,  1972). 

The  Struve  Model 

Origin 

In  1931  Struve  suggested  that  the  emission-absorption 
features  that  are  the  definitive  characteristic  of  Be  stars 
are  formed  in  gaseous  circum-equatorial  rings  surrounding 
(and  ejected  from)  B-type  central  stars  that  are  rotating 
at  their  break-up  velocites.   Later,  Struve  (1942)  was  able 
to  present  the  model  on  the  basis  of  more  detailed  spectral 
considerations.   He  also  proposed  a  circular  velocity  law 
for  the  equatorial  disk, 


V.  (w   ) 
j>c  oc7 


(w  /to   ) 
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where  u  /w    is  the  radial  distance  from  the  central  axis 
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to  any  point  in  the  disk  divided  by  the  stellar  radius  and 

V   (u  )  is  the  rotational  velocity  at  the  stellar  equator. 
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Figure  1,  which  draws  its  inspiration  from  Hack  (1970), 

illustrates  the  model ' s  ability  to  explain  the  major  features 

of  Be  spectra.   The  first  case  is  a  Be  star  whose  axis  of 

rotation  is  perpendicular  to  the  line  of  sight;  the  second, 

one  whose  axis  of  rotation  is  parallel  to  the  observer's  line 
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of  sight.   A  schematic  profile  of  the  expected  emission  at 
one  of  the  Balmer  lines  is  also  shown  for  each  case.   The 
imposition  of  a  small  expansion  velocity  on  the  large 
circular  velocity  gives  a  V/R  asymmetry.  That  the  disk  is 
indeed  drcum  -equatorial  was  given  credence  by  the  work  of 
the  Burbidges.   They  analyzed  (Burbidge  and  Burbidge,  1953) 
high  dispersion  spectra  of  six  stars  whose  emission  features 
were  like  the  second  case  above,  often  called  "pole-on." 
They  found  that,  although  a  very  thin  layer  of  gas  did 
appear  to  exist  above  the  polar  regions  of  the  photosphere, 
the  bulk  of  the  gas  was  in  or  near  the  equatorial  plane. 
Another  aspect  of  the  Struve  model  is  that  it  explains  the 
differences  in  Be  and  shell  spectra  as  due  to  differences 
in  the  density  and  extent  of  the  envelope. 

Although  the  Struve  model  provides  a  qualitative 
explanation  for  the  spectral  features  found  in  Be  stars 
and  shell  stars,  detailed  calculations  of  line  profiles 

require 

(1)  a  model  for  the  radiation  field  of  the  central 

star, 

(2)  a  model  which  gives  the  physical  state, 
density,  and  velocity  of  the  gas  at  all 
points  in  the  disk, 

(3)  the  changes  in  the  radiation  field  resulting 
from  its  passage  through  the  gaseous  disk. 

The  work  of  Mihalas  (1965)  and  of  others  now 
provides  the  first  requirement.   The  second  and  third 
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requirements,  have  not  been  so  satis factorally  achieved  and 
are  hot  even  truly  separable. 

The  Struve  hypothesis  says  nothing  about  the 
mechanism (s)  that  cause  the  formation  and  destruction  of 
the  disk — that  is,  changes  from  B  to  Be  to  shell  star  and 
vice  versa — nor  does  it  say  anything  about  the  mechanisms 
that  produce  other  spectral  variations. 

Line  Profiles 

The  hydrogen-emission  spectrum  is  a  recombination 
spectrum  formed  in  a  disk  the  various  portions  of  which 
Can  move  relative  to  one  another.   A  great  deal  of  work 
on  hydrogen  emission  in  a  stationary  atmosphere  has  been 
done  (Kogure,  1959a,  1959b,  1961).   Sobolev  (1960)  and 
Rothenberg  (1952)  have  studied  the  formation  of  the  hydrogen 
lines  in  spherically  symmetric  envelopes  expanding  with 
constant  velocity.   Marlborough  (Marlborough,  1969  and  1970; 
Marlborough  and  Roy,  1971)   has  calculated  Ha  profiles  using 
a  6-level  hydrogen  atom.   He  used  a  model  for  the  distribu- 
tion and  motion  of  the  material  due  to  Limber  (1964)  which 
will  be  discussed  below. 

Hydrodynamical  Models  for  the  Flow 

Determination  of  the  density  and  velocity  field 
throughout  the  disk  is  a  hydrodynamical  problem.   In  the 
first  such  study  (Limber,  1964) ,a  steady-state  axially- 
symmetric  disk  was  examined  for  a  parameterized  form  of 
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circular  velocity,  law  and  in  the  absence  of  any  radial 
velocity.   Limber  treated  isothermal  envelopes,  poly tropes, 
and  envelopes  which  have  a  specified  temperature  law.   For 
each  case  he  was  able  to  calculate  the  density  at  any 
point  in  the  envelope.   Later,  Limber  (1967)  studied  the 
more  general  case  in  which  a  radial  velocity  component  was 
present. 

Spectral  Variations  and  Their  Possible  Causes 

A  rough  division  of  spectral  changes  in  Be  stars 
will  be  made  here  to  facilitate  the  discussion;  the  divisions 
are: 

(1)  long  term  changes  (time  scale  longer  than 
10  years)  which  characterize  the  appearance 
and  disappearance  of  shell  or  emission  spectra 
as  in  Pleione, 

(2)  medium  term  changes  (1  to  10  years)  particularly 
V/R  changes  with  time  scales  in  this  time 
range, 

(3)  short  term  changes  (  a  year  or  less) . 

There  is  a  certain  inevitability  to  the  long  term 
changes;  Crampin  and  Hoyle  (1960)  argued  that  any  initial 
magnetic  fields  present  at  the  formation  of  the  disk  would 
render  it  unstable  in  100  years  or  less.   Limber  has 
suggested  that  there  may  be  a  non-uniform  transfer  of 
momentum  from  rapidly  rotating  inner  regions  of  the  central 
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star  out  to  the  photo spheric  regions  resulting  in  turn  in 
non-uniform  mass  loss.   Limber  (1969)  analyzed  the  shell 
phase  of  Pleione  in  terms  of  a  time-dependent  mass  flux 
from  the  central  star's  equatorial  regions  to  the  disk.   He 
found  good  agreement  with  the  observations  for  a  flux  whose 
time  dependence  showed  a  slow  increase  to  maximum  followed 
by  a  steep  drop  to  zero. 

McLaughlin  (1961)  reviewed  what  are  called  here 
medium  term  variations.   He  observed  that  the  most  difficult 
problem  concerned  V/R  changes.   McLaughlin  (1962)  discussed 
the  major  attempts  to  explain  the  quasi-periodic  V/R  changes 
of  the  sort  seen  in  ir  Aqr,  and  he  concluded  that  only  a 
suggestion  credited  to  Otto  Struve  was  consistent  with  the 
observations.   Struve' s  suggestion  had  been  that  the  ring 
or  disk  was  elliptical  in  shape  and  that  the  V/R  variations 
result  from  a  line  of  apsides  rotation  of  the  disk.   There 
has  been  no  recent  examination  of  this  idea. 

The  most  comprehensive  work  on  short  term  variations 
is  Huang's  (1972a).   Huang  showed  that  an  asymmetry  in  the 
disk  will  produce  spectral  variations  observable  after  time 
spans  as  short  as  a  fraction  of  a  day;  these  variations  can 
persist  for  as  much  as  a  year  before  the  disk's  differential 
rotation  destroys  the  asymmetry.   Hutchings  (1969,  197Q.) 
suggested  that  photospheric  fluctuations  or  "dumpiness" 
in  the  disk  may  be  responsible  for  short  time  changes. 

Modifications  of  the  Struve  Model 

It  now  appears  that  several  aspects  of  Struve' s 
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original  model  must  be  modified.   First,  Struve' s.  model 
requires  that  the  equatorial  regions  of  the  photosphere 
be  rotating  at  the  break-up  velocity.   Observations  suggest 
that  Be  stars '  rotational  velocities  lie  below  the  break-up 
limit.   Huang  (1972b)  recently  argued  that  the  ejection  of 
matter  is  due  to  both  rapid  rotation  and  a  "temperature- 
dependent  instability."  Huang  feels  that  the  "temperature- 
dependent  instability"  may  be  the  mechanism  studied  by  Lucy 
and  Solomon  (1970)  in  their  work  on  stellar  winds  in  early- 
type  giant  atmospheres;  this  work  argued  for  a  mass  outflow 
due  to  radiation  pressure  from  mid-UV  resonances  in  silicon, 
carbon,  nitrogen,  and  sulfur,  these  being  in  high  ionization 
states. 

Limber  and  Marlborough  (1968)  as  part  of  a  general 
discussion  of  the  physical  processes  at  work  in  Be  stellar 
envelopes,  reanalyzed  the  work  of  Struve  and  Wurm  (1938)  in 
order  to  show  that  the  velocity  law, 
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was  in  better  accord  with  the  data  than  the  circular 
velocity  law  originally  proposed  by  Struve. 

Major  Topics  of  this  Work 
A  principal  contention  of  this  dissertation  is  that 
the  existing  steady-state,  axially-symmetric  hydrodynamical 
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solutions  explain  the  key.  elements  of  the  Be  star  phenomenon. 
More  exact  and  comprehensive  explanation  of  the  body  of 
observational  material  will  require  solutions  to  the  hydro- 
dynamical  equations  which  possess  time  and  angle  dependence. 
The  approach  here  will  be  to  view  the  expression  for  each 
flow  variable  as  the  superposition  of  a  detailed  but  in  some 
sense  second -order  term  onto  the  time-independent  axially- 
symmetric  solution. 

Chapter  II,  a  literature  review,  discusses  the 
previous  hydrodynamical  work  on  Be  stars.   Chapter  III 
studies  the  extent  to  which  the  motion  and  distribution  of 
circumstellar  matter  can  be  treated  as  a  hydrodynamical 
problem.   Chapter  IV  introduces  the  dichotomous  view  of  the 
flow  discussed  above  into  the  hydrodynamical  equations; 
this  chapter  then  presents  the  mathematical  reduction  of 
the  original  set  of  equations  for  the  time  and  angle 
dependence  to  a  set  of  solvable  equations.   The  next  chapter 
(V)  discusses  calculations  based  on  the  solution  to  these 
equations.   In  Chapter  VI  the  energy  equation  is  considered, 
and  Chapter  VII,  after  a  discussion  of  the  relation  of  the 
results  of  this  work  to  turbulence  and  stability  arguments, 
summarizes  the  work. 


CHAPTER  II 
HYDRODYNAMICAL  APPROACHES  AND  STEADY-STATE  SOLUTIONS 

Background 

The  steady-state  axially-syiranetric  solutions  of  the 
hydrodynamical  equations  offer  an  acceptable  description  for 
the  overall  phenomena  presented  by  Be  stars  and  shell  stars. 
These  published  solutions  are  the  starting  point  for  the  work 
which  is  the  principal  subject  of  this  dissertation.   For 
this  reason,  and  because  the  discussion  will  outline  the 
hydrodynamical  approach  to  the  motion  of  circumstellar  material 
in  Be  stars,  the  literature  on  steady-state  axially-symmetric 
solutions  will  be  reviewed  in  some  detail.  After  a  brief 
discussion  of  notation,  the  derivation  and  properties  of 
static  envelopes  will  be  discussed;  then  the  permissible 
solutions  with  non-zero  radial  velocity  will  be  considered; 
finally  the  evidence  available  for  choosing  among  the 
various  possible  circular  velocity  laws  will  reviewed.   There 
will  be  a  short  summary  at  the  last.   The  order  of  presenta- 
tion of  the  three  main  sections  represents  the  order  in  which 
these  ideas  were  developed  by  Limber  and  Marlborough  (Limber, 
1964,  1967;  Limber  and  Marlborough,  1968) . 

The  Hydrodynamical  Equations  and  Notation 
The  hydrodynamical  equations  used  here  are  the  simplest 
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ones  that  could  be  chosen;  they  are  Euler's  equation  and  the 
continuity  equation.   In  dimensioned  coordinates  the  former 
is 
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and  the  latter  takes  the  form 

3p 


__£  +  v  •  (p  v  )  =  0  . 
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The  justification  of  the  use  of  Euler's  equation  will 
be  a  principal  concern  of  the  next  chapter. 

The  coordinate  system  used  here  is  shown  in  Figure 
2.  The  three  coordinates  (to,  $,  z)  may  appear  with  or 
without  the  subscript  c.   The  presence  of  the  subscript 
signifies  that  the  coordinate  in  question  is  dimensioned  in 
cgs  units;  without  the  subscript,  it  is  a  dimensionless 
variable  in  the  system  of  dimensionless  coordinates  which  is 
developed  in  the  next  chapter.  A  similar  procedure  is 
followed  for  the  flow  variables  and  time:  velocity  v, 
density  p,  pressure  P,  and  time  t.  The  equations  already 
presented  in  the  first  chapter  are  consistent  with  this 
practice.   In  cases  where  non-cgs  units  are  used,  these  units 
will  be  explicitly  stated. 

In  the  body  of  this  work  symbols  will  on  occasion 
undergo  some  change  in  definition;  for  instance,  a  variable 
symbol  may  be  required  to  absorb  an  integrating  factor.  This 
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Figure  2   The  cylindrical  coordinate  system  in 
use  here. 
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is  done  only  when  necessary  to  prevent  an  undue  proliferation 
of  closely  related  variable  names.  When  questions  arise 
about  the  meaning  of  a  symbol,  the  best  recourse  is  a  check 
of  t^ie  Key  to  Symbols.  Each  change  in  usage  is  noted  there 
as  is  the  first  page  on  which  the  new  usage  occurs; 

Static  Solutions 
Limber  (1964)  attempted  the  first  hydrodynamical 
study  of  the  material  outside  rapidly  rotating  stars.  He 
took  the  central  star  to  be  rotating  at  the  break-up 
velocity  and  made  the  following  assumptions: 

1)  viscous  and  magnetic  terms  in  the  equation  of 
motion  can  be  ignored; 

2)  radiation  forces  are  either  negligible  or 
includable  through  the  use  of  an  appropriately 
reduced  stellar  mass  in  the  gravitational 
potential; 

3)  disk  self-gravitational  effects  are  negligible; 

4)  steady-state  conditions  prevail; 

5)  only  axially-symmetric  flows  are  included; 

6)  the  Z  component  of  the  fluid  velocity  (perpendicular 
to  the  equatorial  plane)  is  zero; 

7)  the  radial  component  of  the  fluid  velocity  is 
zero; 

8)  the  circular  velocity  has  the  parameterized  form 

v    =  VM°' 

(be      a 

Y        0) 
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where  the  new  symbol  a  is  a  parameter  whose 
value  lies  between  1/2  and  1; 

9)   the  density  is  an  explicit  function  of  the 
pressure  alone. 

The  first  three  conditions  determine  the  nature 
of  the  equation  of  motion  which  can  be  used — Euler's  equation. 
The  third  approximation,  neglect  of  self -gravitational  terms, 
is  quite  reasonable.   The  validity  of  the  first  and  second 
conditions  will  be  discussed  in  later  sections.   The  use  of 
Euler's  equation  as  the  equation  of  motion  is  common  to  all 
hydrodynamical  studies  of  Be  star  envelopes.   The  justification 
for  the  use  of  Euler's  equation  or,  indeed,  any  hydrodynamical 
equation  under  disk  physical  conditions  will  be  the  subject 
of  the  next  chapter. 

The  next  group  of  approximations,  (4)  through  (7), 
reduces  the  equation  of  motion  to  the  form 

2 


vp  -  p.  {  v  (  i  )  +  ^2  I  } 


r_  '     wc  u 


where  va   is  the  circular  velocity  (<J>  component  of  the  veloc- 

<P 

itv)  and  t  is  the  unit  vector  for  the  radial  direction  in 

-*       co 

cylindrical  coordinates.   Conditions  (4)  through  (7)  also 
result  in  solutions  which  wiil  identically  satisfy  the 
continuity  equation,  the  second  of  the  three  equations 
needed  to  specify  completely  the  problem.   Since  violation 
of  (5)  will  show  up  observational ly  as  variations  with  time, 
both  (4)  and  (5)  are  good  approximations  to  the  extent  that 


20 


they  describe  phenomena  with  either  no  time  dependence  or 
long  term  time  dependence.  Near  the  equatorial  plane 
condition   (6)  is  a  good  approximation,  while  condition  (7) 
is  clearly  an  idealization  which  will  be  partially  removed 
in  subsequent  work. 

The  radial  dependence  of  the  circular  velocity 
law  can  be  at  least  roughly  determined  from  dilution  factors 
for  shell  absorption  lines,  from  the  widths  of  absorption 
lines,  and  from  the  widths  of  emission  features.   It  is, 
then,  possible  to  choose  a  parameterized  form  on  the  basis 
of  the  observations.   Calculations  were  done  for  three  values 
of  a  :  ,  1/2,  3/4,  and  1.   In  all  cases  the  circular  velocities 
were  required  to  match  the  stellar  rotational  velocity  at  the 

stellar  surface. 

The  last  condition,  that  the  density  is  a  function 
of  the  pressure  alone,  circumvents  the  failure  to  treat 
the  last  of  the  three  hydrodynamical  equations,  the  equation 
for  energy  transport,  in  the  disk.  The  condition  imposed 
may  be  alternately  expressed  as  requiring  that  surfaces  of 
constant  pressure  and  of  constant  density  be  identical. 

Limber  is  then  able  to  calculate  the  density 
throughout  the  disk  for  isothermal  envelopes,  for  polytropes, 
and  for  envelopes  in  which  the  temperature  though  constant 
along  any  one  surface  of  constant  pressure  may  vary  from 
surface  to  surface. 

The  More  General  Solutions 
The  static  solutions— solutions  in  which  the  radial 
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velocity  is  zero— are,  strictly  speaking,  nonphysical  for  the 
same  reason  that  static  solutions  are  nonphysical  for  the 
solar  wind  (Parker,  1958).  Limber  extended  his  analysis  to 
non-static  solutions  (Limber,  1967) .   Three  assumptions  of 
the  approach  presented  in  the  preceding  section  are  modified; 
they  are  (6) ,  (7) ,  and  (9) .   These  three  are  replaced  by  the 
following  less  restrictive  conditions: 

1)  only  flow  near  the  plane  is  considered; 

2)  the  radial  cross  section  of  flow  must  be  specified 
(in  parameterized  form) ; 

3)  the  temperature  distribution  throughout  the  disk 

is  known. 
As  in  the  static  case  the  ideal  gas  law  is  used.   Euler's 
equation  takes  the  form 


vwc  To£    coc  "   s  3o)c  ^  «c  J   Pc  ^T 


where  the  new  symbol  v^  is  the  radial  component  of  the 
velocity.   The  equation  of  continuity  appears  in  the  inte- 
grated form 

vw(ft»)p(u)  A(w)  =  vu(u0)p(w0)  A(wo) 

where  A(w)  is  the  radial  component  of  the  cross -sectional 
flow.   This  form  of  the  continuity  equation  is  similar  to 
that  seen  in  solar  wind  theory.  The  quantity  A(w)/A(u)0)  is 
analogous  to  the  term 
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in  the  elementary  solar  wind  theory  (Parker,  1963)  where  b  is 
a  parameter  (usually  given  a  value  near  2) .  Limber  chose, 
instead,  a  one-parameter  family  of  cross-sectional  flow 
terms 
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where  p  is  a  parameter.   The  meridional  projections  of  the 
stream  lines  are  straight  lines  running  from  a  point 


0)   =  pw  _ 
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in  the  equatorial  plane  and  sloping  slowly  away  from  the 

equatorial  plane. 

Limber  went  on  to  establish  the  types  of  solutions 
that  are  now  permissible.   Valid  solutions  must  be  subsonic 
below  the  star's  surface  and  approach  interstellar  values 
at  great  distances  from  the  star.   The  only  solutions  which 
satisfy  these  two  boundary  conditions  are  those  whose  radial 

velocity 

1)  first  decreases  with  increasing  <*>, 

2)  reaches  a  minimum  as  w  continues  to  grow, 

3)  then  increases  as  u  continues  its  increase, 

4)  becomes  supersonic  at  sufficient  distance  from 
the  central  star. 

Limber  examined  the  range  of  validity  of  his  earlier 
work,  the  calculations  for  static  envelopes,  in  the  context 
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of  his  new  results.  He  found  that  the  density  calculations 
from  the  earlier  approach  agreed  quite  well  with  the  new 
results  over  most  of  the  near-equatorial  regions  of  the  disk. 
At  large  distances  from  the  star  the  old  solutions  no  longer 
agreed  with  the  new  calculations;  the  distance  at  which  the 
breakdown  of  the  old  approach  occurred  depended  on  the  exact 
physical  conditions/  and  geometric  dimensions  chosen  for 
the  disk. 

Form  of  the  Circular  Velocity  Law 
The  solutions  developed  by  Limber  in  his  19  67  study 
contain  two  parameterized  functions,  the  cross-sectional  flow 
expression  and  the  circular  velocity.   Of  these  two  the  circular 
velocity  is  the  more  crucial;  the  solutions  are  fairly  insen- 
sitive to  changes  in  p,  the  parameter  which  appears  in  the 
cross-sectional  term.   Limber  and  Marlborough  (1968)  examined 
the  physical  conditions  and  observational  evidence  for  the 
circular  velocity's  dependence  on  w  throughout  the  disk. 

The  equatorial  break-up  velocity  is  only  1//2  times 
the  escape  velocity.   This  means  that  if  a  small  element  of 
matter  were  perturbed  with  a  small  radial  velocity  so  that 
it  began  to  move  outward,  and  if  the  element  were  not  further 
disturbed,  then  the  element  would  eventually  fall  back  onto 
the  star.  The  available  evidence  indicates  that  shell  material, 
rather  than  falling  back  onto  the  star,  is  continually  escaping 
the  system;  there  is  mass  loss.   Limber  and  Marlborough  call 
this  the  problem  of  "support."  The  physical  processes  which 
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act  to  provide  the  "support"  profoundly  influence  the  radial 
dependence  of  the  circular  velocity  law.  Limber  and  Marlborough 
divide  the  mechanisms  for  support  into  two  classes:  "direct" 
and  "centrifugal."  The  first  type,  "direct"  mechanisms,  are 
radially  directed  forces  Which  act  on  any  fluid  element  to 
overcome  the  difference  between  the  gravitational  force  on  that 
element  and  the  centrifugal  force.   These  direct  forces  would 
appear  as  additional  terms  in  the  equation  of  motion.  The 
"centrifugal"  classification  refers  to  mechanisms  which  act 
to  transfer  angular  momentum  out  into  the  disk  so  that  the 
centrifugal  force  always  balances  the  gravitational  force  due 
to  the  central  star. 

Direct  mechanisms  do  not  act  to  affect  the  angular 
momentum  of  a  small  fluid  element  which  is  moving  from  the 
stellar  equator  out  into  the  disk;  the  angular  momentum  of 
the  fluid  element  is  conserved.   In  this  case 


V  i=U)„ 

<pc  c 


Quasi-static  motion  of  a  small  fluid  element  away  from  the 
star  with  centrifugal  forces  present  requires  that 


VA  GMo 

to-        2 
c      w_ 


at  any  point  in  the  disk;  hence 


1 
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There  are  four  direct-support  mechanisms  available; 

they  are: 

1)  thermal  support, 

2)  turbulent  support, 

3)  radiative  support, 

4)  magnetic  support. 

Limber  was  able  to  show  that  the  first  possibility,  thermal 
support,  requires  temperatures  and  densities  throughout  the 
disk  that  are  totally  at  variance  with  present  knowledge 
of  the  values  of  these  quantities.  Turbulent  support  requires 
a  highly  turbulent  flow  whose  effects  on  spectral  features 
should  be  observable  but  have  not  been  reported. 

Limber  and  Marlborough  considered  two  possible 
radiative  support  mechanisms,  electron  scattering  and  photo- 
ionization  of  neutral  hydrogen.  The  net  outward-directed 
radiation  force  on  a  small  fluid  element  due  to  Thomson 
scattering  was  calculated  for  typical  disk  electron  densities 
and  was  found  to  be  negligible  in  comparison  to  the  gravitational 
force  acting  on  that  element.   They  noted  that  for  a  given 
neutral  hydrogen  density  there  was  an  upper  limit  to  the  radial 
force  that  could  be  generated  by  photo- ionization  of  neutral 
hydrogen.  Maximum  momentum  transfer  in  the  radial  direction 
from  photo-ionization  occurs  if  the  recombination  following 
the  photo-ionization  emits  a  photon  backwards  toward  the 
source  of  the  original  photon.   To  establish  an  upper  limit, 
Limber  and  Marlborough  assumed  that  all  the  flux  from  the 
central  star  shortward  of  Lyman  a  and  incident  on  the  disk 
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was  absorbed  and  back-emitted.   On  the  basis  of  both  this 
assumption  and  other  considerations,  these  authors  calculated 
an  upper  limit  for  support  due  to  photo-ionization.  This  value 
was  far  less  than  the  value  needed  to  provide  radiative 
support.  While  both  these  calculations  are  dependent  on 
estimates  of  typical  disk  values  for  the  physical  parameters, 
they  are  so  much  smaller  than  the  value  needed  for  direct 
support  that  they  may  be  omitted  from  further  discussion 
until  such  time  as  the  best  estimates  of  these  disk  quantities 
change  dramatically. 

Finally  radial  support  from  magnetic  forces  was 
considered.   The  same  investigators  used  elementary  considera- 
tions to  show  that  the  field  strength,  H,  required  was 

H  >  75  gauss. 

They  then  showed  that  fields  of  this  sort  in  the  disk 

would  disrupt  it  on  a  time  scale  of  days.   The  observational 

results  of. the  disruption  of  the  disk,  loss  of  emission  lines 

or  shell  absorption  lines,  show  much  longer  time  scales. 

Thus  they  concluded  that  magnetic  fields  of  the  required  order 

are  not  present  in  the  disk. 

Analogously  there  are  four  mechanisms  which  can 
act  to  transfer  angular  momentum: 

1)  thermal  viscosity, 

2)  turbulent  viscosity, 

3)  radiative  viscosity, 

4)  magnetic  viscosity. 
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Limber  and  Marlborough  used  the  term  "viscous"  to  denote 
any  phenomenon  which  acts  to  transfer  angular  momentum.   They 
found  that  two  of  these  possibilities,  the  first  and  third, 
were  insignificant  under  disk  conditions.   However,  either 
small  scale  turbulence  or  small  (<  5  gauss)  magnetic  fields 
could  provide  the  needed  angular  momentum  transfer.   They 
concluded  that  the  observational  evidence  was  consistent 
with  either  interpretation.   It  is  important  here  to  note 
that  the  amount  of  turbulence  or  the  magnetic  field  strength 
required  for  angular  momentum  transfer  is  far  too  small  to 
provide  any  direct  support. 

Limber  and  Marlborough  noted  that  the  work  of  Hynek 
and  Struve  (1938)  represented  the  only  attempt  to  draw 
quantitative  inferences  from  observations  about  the  form  of 
the  circular  velocity  law  and  that  this  attempt  had  considered 
only  the  1/w  behavior.   The  data  were  reanalyzed  to  see  if 
the  l//w  form  gave  a  more  consistent  interpretation  of  the 
data.   The  conclusion  of  this  reanalysis  of  the  data  was 
that  the  form  of  the  circular  velocity  law  associated  with 
centrifugal  support  was  in  better  accord  with  the  data  than 
the  old  Struve  form. 

Final  Comments 
Steady-state  solutions  have  been  applied  to  the 
shell  phase  of  Pleione  (Limber,  1969)  by  allowing  the 
envelope  to  be  at  any  moment  very  near  steady-state.  This 
will,  of  course,  work  only  for  long  term  phenomena.  Limber 
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suggested  that  the  end  of  the  shell  phase  represented  the 
result  of  increased  matter  outflow  not  matched  by  sufficient 
energy  flow  into  the  disk;  hence, a  collapse. 

Two  results  of  Limber's  studies  are  of  sufficient 
importance  to  the  remainder  of  this  work  to  bear  restating. 
First,  over  most  of  the  disk,  static  isothermal  envelope 
calculations  agree  well  with  the  more  exact  theory  including 
a  radial  velocity.   Second,  the  circular  velocity  law  is  that 
for  centrifugal  effects.   The  static  isothermal  envelopes 
with  a  circular  velocity  law  can  be  presented  in  closed  form 
and  do  not  require  the  specification  of  a  cross-sectional 
flow  parameter. 


CHAPTER  III 

APPLICATION  OF  HYDRODYNAMICS  TO  THE 
ENVELOPES  OF  Be  STARS 


Disk  Dimensions 


Table  1  gives  the  values  used  here  for  the  physical 
and  geometrical  quantities  which  characterize  the  star  and 
its  envelope.   The  mass  (Mg)  and  surface  temperature  (Tg) 
of  the  star  are  those  of  BO  dwarf.   The  radius  is  larger 
than  one  might  expect  for  such  a  star  and  represents  acknowl- 
edging the  evidence  that  the  star  is  somewhat  evolved  and 
rotationally  distorted.   These  quantities  are  those 
commonly  found  in  the  literature.   The  equatorial  rotational 

velocity  (V  )  represents  a  slight  departure  in  that  the 
s 

rotational  velocity  was  set  at  just  slightly  more  than  half 
the  equatorial  break-up  velocity.   The  disk  temperature, 
which  is  taken  to  be  the  electron  temperature  (Te) ,  the 
mean  molecular  weight  (y)  of  the  disk  gas,  and  the  average 
disk  density  (pJ  are  typical  of  values  found  in  the 
j  ournals . 

From  the  emission  and  absorption  features  at 
hydrogen  and  helium  resonance  lines  one  can  conclude  that 
in  a  typical  Be  star  envelope  a  large  percentage  of  the 
hydrogen  is  ionized  while  most  of  the  helium  is  not. 
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TABLE  I 
Characteristic  Central  Star  and  Disk  Values 


Quantity  Symbol        Value 

Stellar  Mass 

Stellar  Equatorial  Radius 

Surface  Temperature 

Equatorial  Rotational 
Velocity 

Disk  Temperature 

Molecular  Weight 

Average  Disk  Density 


Ms 

10  % 

R 

s 

10RG 

T 

s 

25,000°K 

V 

s 

250  km/s 

T 
e 

10,000°K 

V 

0.68 

P* 

lO"12 
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Theoretical  Foundations  of  Hydrodynamics 
Validity  of  the  Boltzmann  Transport  Equation 

The  most  satisfying  foundation  for  the  hydrodynamical 
equations  is  the  Boltmann  transport  equation.   The  hydro- 
dynamical  equations  are  the  result  of  integrals  of  the 
form 


/  ¥  D  f  dv   . 
1      c  op    C 


*c  is  a  conserved  quantity,  D ,  the  Boltzmann  operator, 
and  f ,  the  distribution  function  which  is  a  solution  of  the 
integro-differential  equation 


op 


Two  assumptions  are  essential  to  the  derivation  of 
the  Boltzmann  transport  equation.   The  first,  the  assumption 
of  molecular  chaos,  is  mentioned  only  for  completeness; 
the  second,  the  restriction  to  binary  forces,  is  open  to 
question  in  the  presence  of  Coulomb  forces.   The  so-called 
collision  integral  of  the  Boltzmann  equation  includes  only 
the  effects  of  binary  collisions.  More  fundamentally,  the 
Boltzmann,  Weltansicht  is  that  the  individual  particles 
comprising  the  gas  spend  most  of  their  time  in  free  flight, 
this  condition  being  interrupted  only  occasionally  and 
briefly  by  collisions. 
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The  long-range  Coulomb  forces  present,  when  charged 
particles  are  constituents  of  the  gas,  generally  result  in 
an  infinite  contribution  from  the  collision  term.   The 
underlying  view  of  the  gas  motion  is  suspect  as  several 
distant  collisions  may  simultaneously  interact  with  one 
another.   There  is,  however,  one  set  of  circumstances  under 
which  the  Boltzmann  approach  is  still  valid  (Zel'dovich 
and  Raizer,  1966) .   The  conditions  are: 

(1)  the  average  Coulomb  potential  energy  at  the 
mean  separation  distance  is  much  less  than 
the  average  thermal  energy; 

(2)  the  Debye  length  is  much  greater  than  the 
mean  separation  distance. 

Assume  that  the  disk  gas  is  made  of  equal  numbers  of 
electrons  and  protons.   The  ratio  of  the  Coulomb  potential 
energy  at  mean  separation  to  average  thermal  energy  for 
such  a  gas  under  disk  physical  conditions  gives 


fff2    -   1.2  xlO"3 

N     kT 
e     e 


where  N  ,  e,  k,  and  Z  are  the  electron  number  density, 
the  electron  charge,  the  Boltzmann  constant,  and  the  average 
ionic  charge.   The  ratio  of  the  Debye  length  to  average 
separation  is 

6.90  (T 0/Nj1/2 

T7T-* *  8>0  ' 
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The  first  condition  is  nicely  satisfied;  the  second,  less 
well  so  but  within  the  limit  of  toleration.   The  conclusion 
to  be-  made  is  that  the  Boltzmann  approach  and,  hence,  the 
hydrodynamical  description  are  justified,  but  are  near  the 
limit  of  validity.   The  additional  species  present  in  a 
more  realistic  representative  disk  gas  will  not  alter  this 
result. 


Relaxation  Times 

A  fundamental  consideration  for  a  gaseous  mixture 
is  the  extent  to  which  that  mixture  can  be  treated  as  a 
single  fluid,  particularly  to  what  extent  a  single  temperature 
can  be  used  to  describe  the  flow.   The  relaxation  time  for 
collisions  between  two  species  is  a  measure  of  the  time 
required  after  an  initial  disturbance  in  one  of  the  species 
to  establish  equilibrium  among  the  translational  degrees 
of  freedom  in  both.   For  a  system  of  the  type  under  study 
here,  there  is  some  time  which  is  characteristic  of  the 
time  scale  on  which  hydrodynamic  variations  occur.   If 
the  ratio  of  the  relaxation  time  between  any  two  major 
species  of  the  gas  and  the  characteristic  time  of  the 
system  is  equal  to  or  greater  than  one,  the  single-fluid 
hydrodynamical  description  is  invalid. 

For  a  gas  under  disk  physical  conditions  there 
are  two  relevant  relaxation  time — that  between  ions  and 
electrons  and  that  between  electrons  and  neutral  atoms. 
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For  purposes  of  studying  the  nature  of  the  disk  gas,  that 

gas  will  be  approximated .by.  a  three-species  gas  containing 

electron,  ionized  hydrogen,  and  neutral  helium;  further, 

the  characteristic  time   (t  )  for  the  system  will  be  taken 

c 

to  be  2irR  /V      . 

The  ratio  of  the  electron-proton  relaxation  time 

(t   )  to  the  characteristic  time  (t  )  is 
e-p  c 


t       3.5  x  10+8V  Tf/2 
e-p  .   s  e 

Tc       2ttR  N  Jin  A, 
s  e   a 


where  A,  is  the  reduced  Debye  length  (Spitzer,  1962) .   The 

ratio  of  the  electron-helium  relaxation  time  (t  .  )    to 

e— ne 

the  characteristic  time  is  approximately 


e-he  s 


2*RsNhe  Te/2  Qhe 


(Zel'doyich  and  Raizer,  1966)  where  Qhe  is  the  helium- 
electron  collision  cross-section  for  typical  disk  conditions 

+4 
and  N,   is  the  helium  number  density.   If  Te  is  10    °K 

—16   2 
Qh  is  approximately  5.7  x  10   cm  . 

For  the  typical  physical  state  of  the  disk,  the 

ratio  of  the  electron-proton  relaxation  time  (xe  )  to  the 
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-4 
characteristic  time  is  4.6  x  10   .   The  ratio  which  compares 

electron-helium  relaxation  time  to  the  characteristic  time 

is  approximately  1.79  x  10""  .   Finally  one  can  compare 


^=^  -  0.4  x  10"1 
Vp 


The  first  two  of  these  three  numbers  indicate  that  the 

temperature  should  be  the  same  for  all  three  components; 

electrons,  protons,  and  neutral  helium. 

These  two  conclusions  should  hold  in  the  more  complex 

disk  gas.   The  principal  disparity  between  the  three-element 

picture  and  the  actual  fluid  is  the  number  of  neutral 

species.   The  quantity  t  .   is  being  used  as  a  measure 

of  the  relaxation  time  for  electron  neutral  collisions.   The 

relaxation  time  for  collisions  of  electrons  with  all  neutral 

atoms ,  in  the  presence  of  additional  neutral  species  besides 

helium,  can  only  be  smaller  than  that  for  helium  alone. 

Therefore,  the  quantity  t  .  ■  .  overestimates  the  actual 

electron-neutral  atom  relaxation  time.   Substitution  of 

a  smaller  number  for  t  .   in  this  section  will  only 

e-he 

strengthen  the  conclusions. 

The  Hydrodynamical  Prescription 

The  Hydrodynamical  Equations 

The  three  conserved  quantities  used  to  generate  the 
hydrodynamical  equations  are  mass,  momentum,  and  energy;  these 
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give  an  equation  of  continuity,  an  equation  of  motion, 
and  an  energy  equation,  in  that  order.  However,  quantities 
such  as  the  heat  flux  and  the  deformation  tensor  which 
appear  in  the  equations  are  defined  by  integrals  that 
contain  the  distribution  function.   Indeed,  the  quantities 
that  are  identified  as  flow  variables  (the  pressure,  density 
and  fluid  velocity)  are  also  so  defined. 

If  the  transport  equation  is  solved  by  a  successive 
approximation  technique  in  which  the  distribution  function 
is  written  as  a  series  of  terms  generated  by  increasingly 
higher  degrees  of  approximation,  a  set  of  hydrodynamical 
equations  may  be  formed  at  each  stage  of  the  approximation. 
The  "zeroth"  approximation,  in  which  only  the  initial  term 
for  the  distribution  function  appears, is  a  locally  Maxwell- 
Boltmann  distribution;  the  resultant  hydrodynamical 
equations  are  those  for  an  ideal  inviscid  fluid.   The 
equation  of  motion  in  this  case,  Euler's  equation,  has  the 
form 


9Vc  +  1 

— £.  +  v   •  Vv  = =-  VP  -  F 

3t  c     C      p     c    c 
c  c 


where  F  is  the  force  per  unit  mass.   The  continuity 
c 

equation  is 


3pc 

— -   +   V  •   (p  v  )  =  0   . 

3tc 
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The  next  approximation  leads  to  the  hydrodynamical 
equations  for  a  viscous  gas.   In  this  case  the  equation  of 
continuity  is  unchanged  but  a  new  equation  of  motion,  called 
the  Navier-Stokes  equation,  is  formed.   The  Navier-Stokes 
equation  will  be  examined  below  to  ascertain  the  extent  to 
which  Euler's  equation  is  a  good  approximation  for  a  fluid 
under  disk  physical  conditions. 

Navier-Stokes  Equations  and  Similarity  Numbers 

The  Navier-Stokes  equation  in  the  presence  of  a 
gravitational  potential  is 


3v_ 


1   r,_     n   „2£ 


+ 


__£  +  v  •  Vv  = —  Vp  +  -^  Vv 

3tc    vc   vvc      pc  vc        pc     c 


c       s 

Pc  rc 


The  quantities  n  and  ?  are  the  first  and  second  vicosity 
coefficients,  respectively.  The  quantities  r\   and  B,   are  both 
always  positive.  The  second  viscosity  coefficient 
represents  effects  that  occur  at  high  density  or  when 
species  with  slowly  excited  degrees  of  freedom  are  present. 
It  is  included  only  for  completeness.  Even  if  these  effects 
were  present,  as  long  as  n  >  £,  all  the  arguments  presented 
below  are  valid. 
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The  second  term  on  the  left-hand  side  in  the  Navier- 
Stokes  equation  is  called  the  inertia  term;  the  second 
and  third  terms  on  the  right  are  the  viscous  terms.   The 
last  term  contains  gravitational  effects. 

Let  RQ  be  a  dimension  characteristic  of  the 
boundary  surface,  Vq  a  typical  value  of  the  fluid  velocity, 
no  a  rePresentative  density.  These  characteristic  values 
are  chosen  such  that  the  associated  physical  quantities 
(rcf  vc,  pc)  vary  from  a  large  fraction  to  several  times 
the  characteristic  values,  that  is,  each  of  the  values 
represents  the  order  of  magnitude  of  the  quantity  of  which 
it  is  a  characteristic.   Characteristic  values  for  the 
pressure  and  time  which  result  from  these  choices  are 


P  =  n  v2 
o    o  o 


and 


T   =  R  /v 
o    o'    o 


respectively.  Note  that  the  characteristic  numbers  are 
dimensional . 

The  characteristic  numbers  can  be  used  to  set  up 
a  system  of  dimensionless  variables  defined  by  the  following 
relations : 

(1)  rc     =  Ror"   , 

(2)  v   =  v  v   , 

c      o 
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(3)  pc  =  nop  , 

(4)  Pc  =  PoP    , 

(5)  t    -  T  t 

CO 

Note  that  the  subscript  zero  quantities  are  dimensioned 
numbers.  These  relations  can  be  used  to  rewrite  the 
Navier-Stokes  equations  in  dimensionless  coordinates.  The 
result  is  shown  below: 


V   „-►     V2  GM 

R0   3t      Ro  °      P    R2      l   r  J 

o 


nVo   vfv  +  (g  +   n/3)   v  V(V  ■  v) 


n  r2    p      n  r2 

o  o  o  o 


2 
Division  by  V/R  yields 


3v  .  ■»■    „-»■    -VP      GMs  /-In 

££  +  v  •  Vv  =  -^- 1-  V  (-  -  1 

3t  p       v2  l  r  i 

o  o 


n    vfv  +   (g  +  n/3)   V(V  ;  v) 


o  o  o   ^      .  o  o  o         K 

Since  the  characteristic  parameters  represent 
the<  order  of  magnitude  of  their  respective  variables,  the 
order  of  magnitude  of  the  ratio  of  any  two  terms  is  the 
ratio  of  their  coefficients.   This  statement  is,  however, 
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open  to  question  near  the  flow  boundary.  The  flow  variables 
may  change  rapidly  over  short  distances  near  the  boundary 
surface  thereby  causing  terms  involving  their  gradients 
to  be  quite  large.   These  coefficient:  ratios  taken  together 
qualitatively  describe  the  flow.   Historically  they  have 
been  called  similarity  numbers  and  given  the  following 
names: 

1.   Reynolds  Number.   This  dimensionless  number, 
labeled  RN  here,  is  the  inverse  of  the  ratio 
of  the  viscous  terms  to  the  inertial  terms, 


RN    R  V  JI 
o  o  o 


2.  Froude  Number.   This  value,  referred  to 
symbolically  as  FN,  is 


FN  =  °- 


Ro(GM/R2) 


Its  inverse  is  the  ratio  of  the  gravitational 
term  to  the  inertial  term. 
3.  Strouhal  Number.   Called  SN,  it  is  given  by 
the  expression 

Vo 

Ro 
The  inverse  of  this  quantity  is  the  ratio  of  the 

time  derivative  of  the  velocity  to  the  inertial 

term. 
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Importance  of  Viscous  Term  in  the 
Equation  of  Motion 

The  relative  importance  of  the  term  in  the  Navier- 

Stokes  equation  can  be  determined  by  calculating  the 

similarity  numbers.   To  calculate  these  numbers,  the 

viscosity  must  be  estimated  and  the  characteristic  values 

(R~r  V  ,  n  )  chosen. 
of  o   o 

Suppose  the  gas  to  be  a  ternary  mixture  composed 
of  hydrogen  ions,  electrons,  and  neutral  helium.   Because 
the  ratio  of  the  electron  mass  to  that  of  either  of  the 
other  two  species  is  so  small,  the  viscosity  of  the  mixture 
is  essentially  determined  by  the  viscosities  of  the 
hydrogen  ions  and  neutral  helium  (Chapman  and  Cowling,  1970) 
For  hydrogen  ions  under  physical  disk  condition  the 
expression  given  by  Spitzer  (1962)  may  be  used.   It  is 


2.2  x  10~15  T5/2 

£nA, 
d 


which  for  disk  temperature  (T  )  gives 


n  =  2.2  x  10"6 


Chapman  and  Cowling  (1970)  give  an  expression  that  can  be 
used  to  calculate  helium  viscosity. 

8a2  W 
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where  a  is  the  molecular  diameter  and  W  is  a  tabulated 
function  ( ir  is  the  mathematical  number)  .   Under  disk 
conditions  the  two  quantities  a   and  W  are 


a  =  2.7  x  10"8  cm2 


W  =  6.00 


This  gives 

n  =  2.4  x  10~4  • 
The  following  characteristic  values  are  used: 


CD   Ro  =  Rs   ' 
(2)   V  -  V  /2ir  , 

O        BJ 


(3)  no  =  pd  . 

For  viscous  effects  due  to  hydrogen  the  similarity  numbers 

are: 

+12 

(1)  RN  =  1.2  x  10 

(2)  FH  =  8  x  10"3  , 


(3)   SN  =  1 


The  helium  viscosity  results  in  a  Reynolds  number  which 
is  two  orders  of  magnitude  smaller  than  that  for  hydrogen, 
but  this  is  of  little  import.  Nor  would  a  more  accurate 
calculation  change  the  basic  result  that  the  Reynolds 
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number  for  the  disk  fluid  is  very  large.  The  ratio  of  the 
viscous  terms  to  the  inertial  terms  is  so  small  that  the 
latter  completely  dominate  the  former.  Removal  of  the 
viscous  terms  reduces  the  Navier-Stokes  equation  to 


a+  „~    GM 

3v 


£  +  v.v*  =  =H      (-1), 


p    R  V* 
o  o 


Euler's  equation. 


CHAPTER  IV 

MATHEMATICAL  ANALYSIS  OF  THE 
HYDRODYNAMICAL  EQUATIONS 


Preliminary  Remarks 

In  the  work  which  follows  magnetic  effects  will  be 
ignored  and  any  radiation  effects  will  be  assumed  either 
entirely  negligible  or  such  that  they  may  be  included  through 
a  small  adjustment  in  the  gravitational  term.  On  the  basis 
of  the  discussion  in  the  last  chapter,  the  viscous  terms  in 
the  Navier-Stokes  equation  can  be  discarded.  The  momentum 
equation  is 


GM 

s 


3v  ^  £   n£  _  _  VP 


3t  P    R  V  ^ 


v  (-  I ) 


The  coordinates,  time,  and  the  flow  variables  are  dimensionless 
and  related  to  the  normal  dimensioned  quantities  through  the 
relation  of  the  preceding  chapter.   In  the  rest  of  this  work 
the  symbol 

GM     GM  (2ir)2 

—    s   _   s 

^  "  R  V2   ~   RV* 
o  o      s  s 

will  be  used.   The  equation  of  continuity, 

|£  +  V  •  (p  v)  =  0   , 
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provides  an  additional  relation  between  pand  v. 

The  full  hydrodynamical  prescription  requires  the 
inclusion  of  a  third  equation,  that  for  energy  transport. 
No  effort  will  be  made  in  this  chapter  to  introduce  such 
an  equation.  Failure  to  study  the  energy  equation  will 
necessitate,  at  some  point,  the  specification  of  an  additional 
condition  linking  two  of  the  flow  variables.   Introduction  of 
this  new  relation  reduced  the  number  of  unknowns  from  five  to 
four.  Euler's  equation,  which  contains  no  viscous  terms,  is 
not  valid  in  the  boundary  layer  between  the  bulk  of  the 
circumstellar  matter  and  the  stellar  photosphere.  No  analysis 
of  the  flow  in  this  region  has  been  made,  nor  will  any  be 
attempted  here. 


Linearization  of  the  Equations  for  Temporal 
and  Angular  Dependence 

Introduction  of  Time  and  Angle  Dependent  Terms 

As  the  steady-state  axially-symmetric  solutions  to 
the  hydrodynamical  equations  appear  to  explain  much  of  the 
Be  star  phenomenon,  each  flow  variable  will  be  each  divided 
into  two  parts:   a  steady-state  axially-symmetric  part  and 
a  term  containing  both  angular  and  temporal  dependence. 
The  notation  is  shown  below: 


(1)   v  =  vo(io,z)  +  V1((l)/^,Zft) 


(2)   p  =  p0(w,z)  +  p1(oj,<J),z,t) 
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3)   P  =  p  (w,z).  +  P1(w,4»fz,.t) 

The  subscript-zero  identifies  the  steady-state  axially- 
symmetric  terms;  the  subscript-one/  the  terms  containing 
the  temporal  and  angular  dependence.   In  order  to  avoid 
confusion  with  the  dimensioning  parameter  P  used  earlier, 
a  lower  case  p  is  used  with  subscript  in  the  case  of 
pressure. 

These  pairs  of  functions  are  substituted  into  the 
equation  of  motion  to  give 


3v    3V.,  ^  ^     + 

"3t  +  IF"  +  Vo  '  Vvo  +  Vl  *  Vvo  +  Vo  *  VV1 


+  Vl  •  Vvx  =  -   (-*-  -  4  +  4  )  VPQ  -  (J:  -  \  )VPl-  V^). 

°    po    po         ^°    po 

and 


3p    3p, 

~£  +  "It   +  ^o  '  Vpo  +  V  Vpo  +  ^o  *  Vpl  +  Vl  *  Vpl  +  poV  *  ^o 


+  P-v  *  vn  +  p,V   •  v  +  p.V   •  v.,  =  0 

O       11        O      1        i. 


The  denominator  in  the  term 


1     V(p„+Pl) 


p0  +  p1      ^o     er 

has  been  approximated  by  first  expanding  in  the  series 
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-1  2    3 

(  1  +  x)    =  1  -  x  -  x  -  x  +  . . 


where  x  is  p,/p  ,      and  then  excluding  terms  in  the  third 
and  in  higher  powers  of  x.   Subsequent  conditions  placed  on 
the  flow  will  limit  the  problem  to  cases  where 

p.  +  p   >  0  . 
Hl   Ko 

With  the  exception  of  this  approximation  for  the  denominator, 
these  equations  are  exact. 

Since  the  subscript-zero  terms  are  themselves  steady- 
state  axially-symraetric  solutions  of  the  problem,  they  satisfy 
the  equations 


r?  +  ;o-  1%°--T0  7po"  V  [~*] 


and 


^|  +  V  •  (P0v0)  =  0 


The  equations  for  the  subscript-one  terms  become 


3v,   +        „        ^         ,Pn   P 

+  v  *  Vv,  +  v,  •  W^  +  v 


o        ro 


°   po 
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and 


ijj|  +  pQV  •  vx  +  PlV  '  VQ  +  Pl  •  Vvx  +  VQ  '  VPl  +  vx  •  Vpo 


+  vx  •  Vp1  =0 


The  underlying  physical  picture  presented  thus  far 
has  been  that  the  steady-state  axially-symmetric  solutions, 
the  subscript-zero  terms,  represent  the  gross  nature  of  the 
flow.   The  subscript-one  terms  are  taken  to  represent  a 
refinement  to  this  gross  nature.   These  refinement  terms 
would  then  be  expected  to  be  small— though  not  negligible. 
The  assumption  will  be  made  here  that  subscript-one  terms  are 
sufficiently  smaller  than  the  subscript- zero  quantities  that 
elements  of  the  equations  which  include  the  products  of 
subscript-one  quantities  can  be  omitted.   The  resulting 
linearized  equations  are 


^o 


and 


!jji  +  pov  •  vx  +  P]Lv  •  vQ  +  S0  •  vPl  +  ^  •  vPo  =  0 
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These  must  be  rewritten  as  scalar  equations 


3t  wl      9oi  wo        9o)  to        3  <J>  u        9<j> 


o«                        3v  ,         2v.    v..             .     ,    9p,        p,    3p, 
qv^  col  4>o   4>1     _        _JL    _    1    .   2i  L 

rzl     ~^T  zo      3z  o)  p„      3oj   •       2      3p 


+     v   ,      ""0      +  v 

3-z     '  fcw      °*  w  Ko      w~        p" 


*VA1  3vA„  3v . ,         v   ,    3v  v        8V, 

3t  oil      3o)  voio      3o)  oi        9<j)  o>        9(f> 


9  3  V0)lV(|)O    ,    VP)QV(j>l 

+     vzl    3z     V(j>o        vzo   3z     V<J)1  o)  o) 


J^_      3?1      +        pl     ^o    , 
po0)      -3?  p2w        ,4 


3v   n  3v_  9v   ,         v.,    3v„^        V.         3v  - 

zl     .  zo  Zl     ,        (pl  ZO     .        (pO  Zl 

3t  Voil        3o)  0)o      3o)  0)        3<J>  oi  3(J> 


3vzo       ..         3vzl    ,       _1  3*1   .    pl  3po 

vzl      3z  zo        3z  p        3z  2      3z      ' 

°  po 


and 


3d-.         d,v  P   v   ,  9p,  3v,%^  9p 

9t  o)  o)  vo)o      3o)  1      3o)  0)1      3oi 


o 


9V    1  1  3pi      •:  3Vd)  3P0  3V<1)1  P13V 

+   po  lS~  +  W      (V<|>o  "9?  +    pl  ~9lT  +  v<j>l  ~9T  +   Po      34>      '    +  ~~3? 


zo 
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•  *PX        3Po  ,.     3Vzl 

vzb  "97  +  vzi  li  +  po 'TT 


Selection  of  a  Subscript-zero  Solution 

The  subscript- zero  solution  which  will  be  used  is 
that  for  a  static  isothermal  envelope  (Limber,  1964) .  This 
solution  is  chosen  for  the  following  reasons: 

(1)  over  most  of  the  disk,  the  isothermal  envelopes 
are  in  agreement  with  calculations  based  on 
more  advanced  treatments  (See  Chapter  II) ; 

(2)  the  solution  is  in  closed  form; 

(3)  the  simplicity  of  this  solution  results,  in 
turn,  in  a  comparatively  simple  set  of  expressions 
for  the  subscript-one  terms. 

This  subscript-zero  solution  (in  dimensionless  coordinates) 
is 

v   =  0   , 

wo 


v   =  0   , 

zo      ' 


1  rl 

P 


=  exp  f  -ff:fe"  J?  ♦  ,i,V*  ]   >' 
Po  -  Vexp  {-  §  (i  -   2  +  ^j  )  I  • 


The.  quantity  Q  is  defined 


kR  T 
n  -   a  e 

u  ~  ymM  G   ' 
s 
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while  f  is  given  by 


.  .(.2.7f.)2kt 


um  Vs 


The  ratio  of  these  two,  Q/*,  is  the  Froude  number  (FN) . 

Substitution  of  these  expressions  into  the  linearized 
equations  for  the  subscript-one  quantities  gives 


SV»1     +       2n      8V"1  -     _  ... 

-sr-  +  ^375-  -sir    ^n  v*i 


-T7I    »«--*!-•[§  1^44     ^ 


VPl  {exp[  |  (  i  -  i  }]}^  {exp[-|  t  i  -  |)  ]  }      . 


#    -  -fr  ^^?-J<J|(:-i)^ 


3vzi   . .  _w  !l5i . .  ,exp  ri  ri .  1 1 1 )  !!i 

TT"  "37?     8*  '     P  LQ   1m       rJJ'        32 


ID 


+  PlMexp[|(i-i)]la|C-i(i     §)]»    - 

and 

!^+ :j5a  {exp[-  i  (|-  i)  ]  }  +  vul  ^  t  exP[  -  §(i-  i)  ]} 
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+  I  '•»*[-  §  (I  -  I  )]J  a!  V  +  "§72  ^  +  =  <expt  "ffCs  "  ?']  » 


^Ujjji     {exp    [-§(i-i)    ]}}vzl 

Mexp[-i(I-i)],4|i=     0      . 

While  the  formal  division  of  the  flow  into  two 
elements  is  quite  general,  the  division  has  physical  significance 
only  if  what  are  called  here  the  subscript- zero  terms  adequately 
represent  the  steady-state  axially-symmetric  part  of  the  flow. 
The  static  isothermal  solutions  do  not  adequately  represent 
such  flow  at  large  distance  from  the  equatorial  plane.  The 
solutions  to  the  above  equations  are  useful  only  for  regions 
near  the  plane.   It  will  be  the  practice  in  the  rest  of  the 
chapter  to  limit  the  discussion  to  the  region  |z|  /w  £  0.1 
(quite  a  conservative  value) .   Such  a  limitation  may  seem  a 
drastic  approximation;  in  fact,  it  is  not.   This  "thin" 
region  is  a  solar  diameter  in  thickness  near  the  stellar 
surface  and  larger  farther  out.   Additionally,  examination 
of  the  z-component  of  Euler's  equation  reveals  that  only  the 
pressure  gradient  balances  the  z-component  of  the  gravitational 
force,  but  the  ratio  of  the  former  to  the  latter  gives 


-1  P 

-2-Z »     i±8.3xl0"3 


^or  (?) 


qo 
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The  dominance  of  the  z-component  of  the  gravitational  term 
suggests  that  the  disk  material  is  highly  concentrated 
toward  the  plane. 


Considerations  Concerning  Approximations  Near  the  Plane 
For  z  <  to  the  term  1/r  becomes 


^-^[i-ci)2]-3/2-^[i-!(!Afc§>4...] 


--Ji  [i-l  f^)2]  - 


similarly,  the  1/r  term  is 


i-I  [l  -  I  f  £)21 
r   u  <-  *   2  l  oi  J   -1 


All  derivatives  of  known  functions  should  be 
performed  before  approximations  such  as  those  described 
above  are  included.   This  has  been  the  practice  here.   On 
substitution  of  the  appropriate  expressions  for  1/r3  and 
1/r,  the  set  of  four  equations  becomes 


JMl  +    2tt   .  ■  oil     4tt        p     t        z2  ^  1  3pl 


3*Pt     2   r        2 
2Q    io*   L      2Qu3   J 
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8>        +   _jr  +  _2*       J^l  _  .   1     f exp    r_zL  )  1    !!l 


3vzl        ,       2ir       8vzl  ^  -   -*     -  ^        *^z    -  -   -2 


TF 


♦  ifr -#-..- [«-&)]*--$  [~.(5y] 


and 


ZVzl       r  r         z2n    J  f       z2     .1   9vzl  _  n 

Qo)  2Qw  2Qoj 


The  form  of  the  last  and  most  complicated  of  the  set 
of  four  equations  suggests  a  change  of  variable  which  will 
simplify  the  equation.   Define  the  variable  v^^  by  the 
expression 

e  r  z   i 

vil  =  Vil  eXp  ( 3  J  ' 

XX     1J-  2QuT 


the  i  standing  for  o),<{>,  or  z.   Now 
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2  9vA,         _  2      „  3v 


♦*[•*  t-ii)]  -#  ♦  c  -  <- ±*  9 


zl 


20V       J         ^  L  2Quf 


3z 


yi  j    exp j  J 

Qui  L  2Q   J       J 


Comparison  of  this  result  with  the  terms  in  the  last  of  the 
set  of  four  equations  above  shows  that  the  last  equation  can 
be  written 

Study  of  the  first  three  members  of  the  set,  which  were 
derived  from  Euler's  equation,  indicates: 

(1)  that  no  component  of  v-j^  appears  in  a  partial 
derivative  with  respect  to  either  u  or  z, 

(2)  that  the  exponentials,  when  they  appear  in  any 
term,  come  before  all  operators. 

On  multiplication  by  exp  (-z2/2Qw3)  and  with  the  introduction 
of  the  variable  v?  ,  the  first  three  equations  of  the  set 
become 

col    -   col  _  9f  e  ..  _  ±   +  _  p. 


■•*»5i..-ftCi  ,   f3vn 

1 

3px 

3t             2                    3d) 

0) 

3(J) 

and 
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-9t~      9*       3z    Qw3 


Note  that  the  function  f  used  above  is 


_   2tt 
f   "372  ' 


The  exponential  in  the  definition  of  v^  makes  the  function 
strongly  peaked  toward  the  equatorial  plane  . 


Additional  Relations 

There  are  four  equations,  but  five  unknowns.  As 
mentioned  in  the  beginning  of  this  chapter  the  missing 
equation  is  that  for  energy  transport.   This  equation 
requires  more  physical  information  about  the  nature  of  the 
flow  than  is  needed  for  either  the  equation  of  motion  or 
the  continuity  equation  and  it  is  mathematically  more 
complex  than  either.   The  traditional  approach  has  been  to 
circumvent  the  need  for  such  an  equation  by  specifying 
in  some  manner  the  relationship  between  any  two  of  the  five 
unknowns;  this  will  also  be  done  here.   The  density  and 
pressure  will  be  assumed  to  be  related  by  a  condition  of 
the  form 

Pi  -  rPn 
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where  T  is  a  parameter.  As  an  illustration,  suppose  the 
relationship  between  the  subscript-one:  density  and  the 
subscript-one  pressure  to  be  isentropic  with  the  isentropic 
constant  determined  from  the  subscript-zero  quantities, 
then 

2.  , 

r  =  — ■£ — 7  .  i<T2 

(vg/2  )2 

where  c  is  the  isothermal  speed  of  sound  determined  from 
the  subscript-zero  terms. 

The  boundary  conditions  for  viscous  hydrodynamics 
require  that  both  the  tangential  and  normal  components  of 
the  velocity  on  either  side  of  the  boundary  between  flow 
and  bounding  surface  be  equal.   In  the  steady-state  axially- 
symmetric  solutions,  the  interior  boundary  surface  is 
rotating  and  has  no  radial  velocity.   The  steady-state 
solutions  do  not  represent  a  true  use  of  boundary  conditions, 
for  the  velocity  dependence  is  in  effect  assumed.  Further, 
the  flow  is  taken  to  be  Eulerian  up  to  the  boundary.  An 
adequate  treatment  of  the  boundary-value  problem  requires: 

(1)  treating  the  flow  in  the  boundary  layer, 

(2)  describing  the  behavior  of  both  the  tangential 
and  normal  components  of  the  velocity 

over  the  boundary  surface,  particularly 

dependence  on  time  and  angle. 
There  is  no  adequate  treatment  of  the  boundary  layer  flow, 
and  the  observations  provide  little  information  about  the 
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temporal  or  angular  dependence  of  the  stellar  surface 
velocity.   For  these  two  reasons  boundary  conditions  are 
of  little  use  in  the  study  of  the  nature  of  the  subscript-one 
solutions . 

To  this  point  the  four  equations  which  together 
comprise  Euler's  equation  and  continuity  are 


8V..,      3v..,  3Pi      3¥.z2 

2Qu> 


1         DV  ,  "HT 

0)1    ,    f  0)1        r,f       „        _        V       ±  +  -. 

3t~  +-  f  ~W~  V'*l     ?   3a)     ^^   : 


!lii  +  f  va>i  +  !!l*i  -  _  i  !!i 

~9t    1 —    3*   "    w   3<J> 


3v,t      3v,t       3PX   VZP]. 


and 


3p,         3p,         v  ,         3v   ,         ,    3v,.,         3v. 
_i+f  —I  +  J^  +        "1  +  I        ft1  +        zl  =  o 
3t    *     3*  u>  3u>  w      a<fr  3z 


Note  that  the  e  superscripts  have  been  dropped.   The  symbol 

->- 

v,  now  contains  the  exponential  term,  and  will  for  the  rest 

of  this  chapter. 

Integration  Over  z 
Further  Restrictions  on  the  Flow 

The  subscript-one  flow  will  be  assumed  to  possess 
symmetry  about  the  equatorial  plane;  that  is,  the  flow  below 
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the  plane  is  the  mirror  image  of  the  flow  above  the  plane. 

This  condition  requires  that  the  flow  variables,  with  the 

exception  of  v  , ,  must  be  odd  functions  of  z . 

The  set  of  four  equations  which  describe  the  flow 

will  be  integrated  with  respect  to  z  over  the  interval 

r-z  ,  +  z  ]  where  z   is  a  function  of  to  alone.   New  variables 
1  o'    o        o 


must  be  introduced;  they  are: 

+z 

p 

-z 


(1)   ux  -  /  °  V;^  , 


+zo 


(2)  a1   =  /    P]dz   . 


-z 
o 


Under  the  condition  requiring  symmetry  about  the  equatorial 
plane,  the  z-component  of  the  integrated  velocity  (ulz)  xs 
identically  zero.   The  integration  with  respect  to  z  and  the 
introduction  of  new  variables  require  that  operator  reversals 
of  the  type 

+zo(w)  f+zo(w) 

'-■o0l)  l£  g(a)'Z)  dZ  ■  TO    -Zo(03)g(a)'Z)  dZ' 

where  g(co, z)  is  a  well-defined  function,  be  performed.   The 
justification  for  such  changes  in  the  order  of  operators 
appears  in  Appendix  A. 

These  two  steps  result  in  the  three  equations 

W  %1+  f  ST  ".1  -  2fV  -  r  TS  +  ^T    L°o  Z*  PldZ     ' 
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3u<j,i    x    f  .3  ■  r  dai 

Tt  ■ .       ?  uwi  +  f  1?  U<*>1  "  »  "aT    ' 


and 


9a,        U .,        .  ,.  -, .  .Q  Bo-, 

1  +  -J2i  +•■-£.  u  ,   +  -  It  u*i   +  f  T3~-  +  2v*i  <z  =  zn>   "  ° 

3t  w  •  •     9a>     wl        to   3c{>      $1  3 9  zl  o 


The  last  term  of  the  first  equation  above  will  be 
rewritten 


/   °  z2p,  dz=A/  °  p.  dz 

-z  -z^ 

o  o 


where  the  quantity  A  may  be  considered  defined  by  this  rela- 
tion.  A  is,  strictly  speaking,  a  function,  but  it  will  be 
treated  as  a  constant  parameter.   The  validity  of  this 
approximation  can  be  tested;  the  solutions  should  change 
only  slowly  with  changes  in  A. 

The  last  term  of  the  fourth  equation,  2vzi(z0)' 
represents  mass  flow  perpendicular  to  the  equatorial  plane; 
a  quantity  generally  believed  small.   It  is  a  reasonable 
approximation  to  set  this  term  to  zero.   Such  a  step,  while 
not  crucial,  does  simplify  the  calculation. 

The  third  term  in  the  continuity  equation,  which 
contains  a  partial  derivative  with  respect  to  to,  should 
be  dominated  by  the  neighboring  terms  which  contain  <f)  and  t 
partial  derivatives.   This  term  will  therefore  be  dropped 
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from  further  consideration, as  is  a  similar  term  in  the  first 
equation  of  the  set. 

The  result  of  all  these  steps  is  shown  below: 


(!)   |J  a,  +  f  {£  m  -  2f  u. .  =  2™      < 

3t   1     3<j>   1  <)>1    2Qto4   1 

/ox  3u   ,  f     _  3u  .  _   r   l 

<2>   3t  *1  +  2  VL  +  f  9*  *1  "  "  S  3T   » 

and 

3a,  3^   ua)  x  tofa 


Solutions 

The  following  symbols  are  introduced  to  simplify  the 
notation : 

(1)   cl   ~~2Q   » 


(2)  c2  =  -r  , 


<3>  L  "  it  +  f  I? 


The  equations  of  the  previous  section  become 


L  uu,l  "  2*  V  "7"   ' 


f        C2;  ?°1 

L  U((>1  +  2  uul  '*  u  "Tf  ' 
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and 


Lff    +  JSk  +  i  __M.  =  o 

1         -a)  0)      3<j> 


This  set  of  simultaneous  differential  equations  must 
now  be  reduced  to  a  set  in  which  there  is  a  separate  equation 
for  each  independent  variable.   This  involves  much  differentia- 
tion and  manipulation.   These  reductions  are  shown  in  Appendix 
B.   The  results  of  these  manipulation  are  that  all  three 
equations  share  the  common  form 

(variable)  =  0 

where  "variable"  may  be  either  u^/w  ,  u^  ,  or  a±    . 
Each  of  these  three  equations  has  the  form 

»  (i&  '  st  •  •>  y  -  o   ■ 

Let 

y  ■  exp  [a(u))<j>  +  S(a>)t] 

where  a(u)  and  3(a))  are  as  yet  unspecified  functions.  Since 


i  exp  {a(w)<(>  +  3(w)t]  «  a(w)  exp  [a(u))<{>  +  3  (w)  t] 
3<J> 
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and 


|r  exp  [0(w)*  +  3(w)t]  -  3  (to)  exp  [a(w)<j>  +  3(to)t], 


one  has  that 


F  (|j  ,  -j|  ,w)   y  =  F(a(a)),  B(w),u)y 


A  solution  (y)  exists  only  for  a  (go)  and  3  (to)  such  that 

F^a(u),  3(u>),w)   =  0. 

Unfortunately  there  are  infinitely  many  such  (a  (to),  3 (w))  pairs. 

The  temporal  and  angular  dependence  of  the  flow 
variables  takes  the  form  exp[ikt  +  in<j>]  where  the  substitutions 

a  =  ik" 

3  =  in 

have  been  made.  The  quantity  ic  is,  in  effect,  a  frequency. 
For  comparison  with  observations  the  associated  quantity 
2Trt/k  ,  which  is  the  period  associated  with  such  a  frequency 
is  more  useful.   Since  $   and  (<j>  +  2tt)  are  physically  the 
same  point,  n  must  be  an  integer.   The  equation  which  k  and 
n  must  satisfy  is  derived  in  Appendix  B.   This  equation  is 

c  n  c 

ok3  +  3wnfk"2  +  ( — wf 2  +  3eofi 2n2  -  -|  )  k" 
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c2n3f   2c  fn     -      3  ,   c^nf 

(D    -    co  •   2w 


For  a  given  n  and  u  there  will  be  in  general,  three  values 
for  k~. 


Relations  Among  the  Flow  Variables 

The  temporal  and  angular  dependence  of  each  flow 
variable  takes  the  form 

exp  [in  (J)  +  ifc~(oj,n)t] 

The  quantity  n  is  an  integer,  and  the  Ict^n)  are  roots  of  a 
polynomial  equation  whose  coefficients  are  functions  of  w 
and  n;  The  general  solution  for  the  problem  may  be  written 


u  ,  =   7  U  i-(u)   exp[in<{>  +  ik(w,n)t], 
wl   n~k  nk 


u*l  =  E-CnJc^   exp[in<j>  +  ik(w,n)t], 
*         n,k 


and 


<*i  m     I  S  7-(w)   exp[in<}>  +  ik(to,n)t] 
1   n,k  nK 


The  relationships  among  the  three  quantities,  U^,  C^  ,  and 
S  r,  determine  the  relationship  among  the  flow  variables. 
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Substitution  of  the  general  solution  into  the 
original  set  of  three  equations  yields 


7   [i(k  +  nf)  Un£  -  2f  Cn5r  \±   S^]{  exp[in<f>  +  ik(o>,n)t]}  =  0 
n;k  4 


J _[  _2*  +  Si  +  i(k  +  nf)  S^  ]{  exp[in<j>  +  ik(aj,n)t]}  =  0/ 

n,k   2       w 


£_[_  -4p  +  i(k  +  nf)  Cn^  -  -^  Sn^J{exp[in<j.  +  ik(u>,n)t]}  =  0- 
n,k 


These  three  equations  all  require  summations  over  n  and  k. 
What  is  desired  is  a  set  of  relations  among  the  three 
quantities  (U  r,  Snp  C  ^)  for  a  given  n  and  k. 

Each  equation  is  multiplied  by  the  quantity 

exp[in'<J>  +  ik'  (w,n')t]  . 

Next,   double  integration  of  the  form 

,  +T"  +TT 

^  Lim    — i_     /  dt      /        d<J>  exp[-in<}>  -     ikt] 

T   -v  oo     4Ut  _ 

-T       -IT 

is  performed  on  each  equation  (and  the  order  of  summation  and 
integration  reversed) .   One  finds 

i(k'  +  n'f)  Un,p  -  2f  Cn.K,  -  '^   SnlJ-,  =  0  , 

0) 
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in'c. 


|  Un,£,  +  1(5  +  n'f)  Cn,p  -  — ^  Sn^,  =  0 


» 


SVP  +  ¥   Cn'k"'  +  i(k'  +n'f)  Vic' 


In  the  following  pages  the  primes  will  be  surpressed . 

These  expressions  form  a  system  of  homogeneous  linear 
equations.  There  is  no  a  priori  reason  to  assume  that  any 
non-trivial  solutions  exist.   If  solutions  do  exist,  they  will 
constitute  at  best  a  single  infinity  of  solutions.   Such  a 
case  can  only  provide  unique  values  for  the  ratios  of  two 
of  the  unknowns  with  respect  to  the  third. 

If  non-trivial  solutions  exist,  the  determinant 
formed  by  the  coefficients  of  the  variables  in  the  system 
above, 

i(k  +  nf)  -<2f  -c^/w4 

f/2  i(k  +  nf)  -in  c2/oj 

1/w  in/u  i  (k  +  nf ) 

should  be  zero.   The  details  of  the  calculation  of  the 
determinant  are  left  to  Appendix  C  but  the  results  are 
pleasantly  familiar.   The  condition  that  the  determinant 
be  zero  is  a  polynomial  in  k  which  must  be  zero.   This 
polynomial  is  the  same  as  that  used  to  calculate  k.   Therefore 
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the  determinant  is  identically  zero. 

There  will  be  a  single  infinity  of  solutions  if 
the  rank  of  the  matrix  from  which  the  determinant  above  was 
formed  is  two;  then  the  equations  can  be  solved  for  the 
ratios  of  two  of  the  unknowns  with  respect  to  the  third. 
The  rank  of  this  matrix  is,  indeed,  two  (See  Appendix  C) . 
The  ratios  which  will  be  determined  are 


and 


Snk 


They  satisfy  the  system  of  equations 

-  (k  +  nf )  x  -  2fy  =  -i 


°1 

or 


fv  mc0 

"  +  i(k  +  nf)y  =  -j-S. 


I +  HF -!«♦»« 


These  two  ratios  are 


re.  (k  +  nf)   2fnc0  -i 

*  -  i|_ 1 +  "  J 


w     .      0) 
DET 
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r  nc2(]c  .+.  nf)   .  c^f  i 


DET 


where 


DET  =  [-(k  +  nf)2  +  f2]  . 


A  Reprise 
The  analysis  of  this  chapter  is  valid  only  to  the 
extent  that  the  underlying  physical  picture,  discussed  in 
the  second  section  of  this  chapter,  is  a  reasonable  representa- 
tion of  Be  stars.  This  requires  the  dichotomous  view  of 
the  flow  presented  there  to  be  a  real  one.   The  steady-state 
axi ally- symmetric  effects  must  dominate  the  flow.   The 
elements  of  the  flow  which  depend  on  time  and  angle  must 

be  secondary. 

Further , solutions  were  obtained  only  when: 

(1)  known  functions  appearing  in  the  differential 
equations  were  approximated  by  formulae  valid 
only  for  z  <  u; 

(2)  the  flow  variables  were  assumed  to  possess 
symmetry  about  the  equatorial  plane; 

(3)  two  parameters,  r  and  A,  were  introduced; 

(4)  the  equations  were  integrated  with  respect 

to  z  (and  new  dependent  variables  introduced) ; 
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(5)  two  terms  in  the  continuity,  equation  and  one 
in  the  equation  of  motion  which  appear  to  be 
of  secondary  importance  were  dropped. 

The  solutions  to  this  set  of  equations  are  proportional 
to  exp[in<J>  +  ik(a),n)t],  where  h  is  any  integer  and  the  values 
of  k  for  any  n  and  w  are  the  roots  of  a  polynomial  equation. 
The  general  solution  for  the  flow  can  be  determined  to 
within  the  infinite  set  of  multiplicative  factors  S  ■«  which 
appear  in  the  integrated  density.   Removal  of  the  S  ■?- 
indeterminance  awaits  a  boundary  value  analysis. 


CHAPTER  V 
.COMPUTATIONS  AND  RESULTS 

Introduction  to  the  Computations 

Both  the  angular  and  the  temporal  dependence  of  the 
solutions  in  the  previous  chapter  are  controlled  by  the 
integers  n  and  the  quantities  k,  respectively.   The  calculation 
for  k",  which  depends  on  n,  u,  A,  and  r,  emerges  as  the 
principal  difficulty.  Accordingly,  the  three  principal 
concerns  of  this  chapter  will  be: 

1)  the  nature  of  the  roots  of  the  equation  for  k, 

2)  the  means  used  to  calculate  them, 

3)  the  results  of  the  computations  for  various 
choices  on  which  the  k  equation  depends. 

As  the  solutions  cannot  be  fully  evaluated  in  the 
absence  of  a  proper  boundary  treatment,  these  results  will 
constitute  much  of  what  can  be  learned  from  the  analysis. 
In  most  of  the  discussion  the  period 

2ttT 
IperI  = 

|k"| 

(generally  in  hours)  will  be  used. 

The  first  major  section  of  this  chapter  is  devoted 
to  the  first  two  topics  above,  the  nature  of  the  roots  and 
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the  procedure  used  to  compute  them.   The  results  of  these 
computations,  the  third  major  concern  of  this  chapter, 
naturally  group  themselves  into  three  divisions: 

(1)  calculation  of  the  pertinent  quantities 
in  the  absence  of  any  angular  dependence 
(n  =  0), 

(2)  temporal  behavior  in  the  presence  of  the 
simplest  angular  dependence  (n  =  1) , 

(3)  the  trends  of  the  results  for  more  complex 
angular  behavior   (n  <  1) . 

The  three  central  sections  of  the  chapter  reflect  this 
three- fold  division.   In  each  the  dependence  of  the 
results  on  <of  A,  and  r   is  discussed.   The  penultimate 
section  presents  sample  calculations  of  the  flow 
variables  which  illustrate  the  qualitative  aspects  of  the 
solutions.   The  last  section  of  the  chapter  compares 
the  results  with  relevant  observations. 


The  Equation  for  k 
The  Consequences  of  imaginary  Roots 

The  equation  for  k,  derived  in  Chapter  IV  is 

2 

uk.3  +   3umfk2  +    (  r^ wf2   +  3ton2f2  -  -J  ")      *" 

0)  to 
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c,n3f   2c  nf     3      ,  3   c^nf 

a)       a)  2co* 


Being  an  equation  of  the  third  degree,  this  equation  has 

either  three  real  roots  or  one  real  and  two  complex  roots. 

If  there  are  complex  roots,  one  is  the  complex  conjugate 

of  the  other.   The  appearance  of  complex  roots  had  immediate 

and  serious  consequences  for  the  methodology  of  Chapter  IV. 

One  of  these  roots  must  of  necessity  have  a  negative  imaginary 

part  which  will  appear  in  the  solutions  as  an  increasing 

exponential  term,  exp  HE^t),  where  Ik-J  is  the  absolute 

value  of  the  imaginary  part.   Flow  variables  which  increase 

exponentially  in  time  will,  at  some  time,  violate  the  condition 

under  which  linearization  of  the  flow  variables  is  valid. 

It  is  possible  to  choose  the  three  multiplicative  factors 

(S  T-,   C  n  U  r)  in  such  manner  that  these  increasing  terms 
v  nk'   nk'   nk 

do  not  contribute  to  the  final  solution,  but  such  choices 
are  suspect.   Only  for  conditions  in  which  all  the  roots 
are  real  can  the  linearization  procedure,  and  hence  all 
thereafter,  be  considered  valid. 

Initial  Values  for  the  Parameters 

The  quantities  r  and  A,  which  appear  in  the 
coefficients  of  the  equation  for  k"  (in  c,  and  c2  respectively), 
enter  the  problem  through  the  relations 

Pi  =  rpl 
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and 


+z  +z 

J     °  zz  P±   dz  =  A  /  °  px  dz   , 

-z  -z 

o  o 


respectively.  The  calculations,  whose  results  are  described 
in  later  sections ,will  be  performed  over  a  whole  range  of 
values  for  each  of  these  quantities;  however,  an  initial 
or  representative  value  for  each  quantity  must  be  determined. 

The  representative  value  for  r — called  T — is  the 
value  for  r  in  the  case  that  p  and  p.,  are  related 
isentropically .   The  isentropic  coefficient  is  determined 
from  the  subscript- zero  quantities;  it  is 


(2ir)2  kT 

r {-^ 

v2   ym 

s 


For  an  electron  temperature  (T  )  of  10   °K, 


r  =  7.67  x  10"2 


A  representative  value  for  A — called  A — can  be 
calculated  be  setting  p1  equal  to  a  constant;  for  z  equal 
to  0.1000, 

A  =  -2  =  3.334  x  10~J   . 
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Machine  Computation  of  the  Roots 

Although  the  mathematics  required  to  determine  the 
values  of  ic  at  any  point  in  the  disk  is  not  difficult,  it 
is  time  consuming.  This  calculation  must  be  performed  at 
many  points  and  for  a  number  of  different  values  for  the 
quantities  n,  A,  and  T;  machine  computation  is  neccessary 
to  perform  as  many  calculations  as  are  needed.  A  program 
was  devised  which  will  perform  the  following: 

1.  Input.   The  program  reads  in  the  values  which 
are  to  be  assumed  by  the  physical  and  geometrical  quantities 
appearing  in  the  coefficients  of  the  ic  equation.   The  program 
must  also  be  given  initial  values  for  A  and  r. 

2.  Roots  of  the  k  equation.  On  the  basis  of  internal 
control  statements,  the  program  establishes  a  sequence  of 

A,  r,  n,  and  u  combinations.  For  each  of  these  combinations, 
the  program  calculates  first  the  coefficient  and  then  the 
roots  of  the  equation  for  k.  Each  root  is  searched  for  an 
imaginary  part;  if  the  root  is  real,  the  quantity  2Tr/k  (in 
hours)  is  calculated. 

3.  Solutions  for  the  flow  variables.  When  given  S  i- 
as  a  function  of  u,  the  program  uses  the  formalism  developed 
in  Chapter  IV  to  calculate  both  C^  and  U^.   The  quantities 


Re{Un^(oj)  exp[in<J>  +  ik(tu,n)  t] } 
Re{Cnk((t))  exPtin*  +  ik(u),n)t]}   , 
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and 

Re{S  -r  (w)  exp[in<t>  +  ik(ai,n)t]} 

are  calculated  by  the  program  throughout  the  disk  for 
various  choices  of  angles  and  times. 

4.  Readout.   The  program  prints  the  results  of  the 
calculations  in  tabular  form. 
The  actual  program  steps  are  contained  in  Appendix  D. 

The  Case  n  =  0 


Distinctive  Aspects  of  the  n  =  0  Case 

The  k  values  here  are  for  solutions  with  no  angular 
dependence.  In  the  absence  of  dependence,  the  k  equation 
takes  the  form 

wk3  -  (uf2  +  -i)  k"  »  0   . 

the  three  roots  are 

o   c  1/2 
k  =  0,  +  K±     +  -3) 

The  first  value  corresponds  to  a  trivial  solution.   The 
last  expression  is  real  for  all  reasonable  values  of  c^. 

Variation  with  coy  r  and  A 

The  expression  for  the  non-trivial  roots  depends  on 
a)  and  c,  but  not  on  c_.  Now 


Figure  3   The  dependence  of  |PER|  on 

to  is  shown  for  the  case  n  =  0 
and  c,  =  36.122. 
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-^-2  =  0.152 


where 


(2ir) 


-  _  3M 
1    2Q 


Setting  c,  =  c.   one  finds  that 


2tt 


(  i+  (~^t)  ^  +  .-) 


|k|  =  "^  i**  ^I^7J? 


U) 


~   2tt 
"  ^72  ' 

a 
The  absolute  value  of  the  corresponding  period  is 


or 


PER|=  27T.   T 
Ikl    ° 


|PER|=  48.5   w3/2 


For  large  c,  the  behavior  of  |per|  near  the  stellar  surface 
(<o  =  1)  departs  slightly  from  w  '   behavior,  as  would  be 

expected  from  the  full  expression  for  |k~|  ,  but  the  outer 

3/2 
regions  of  the  disk  still  show  the  a)  '   behavior.   These 

features  are  shown  in  Figure  3. 


The  Case  n  =  1 
Preliminary  Remarks 

The  n  =  0  case  cannot  be  used  for  a  general  discussion. 
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Exploratory  calculations-indicated  that  the  n  =  1  case 
illustrates  most  of  the  characteristic  features  of  the 
general  non-zero  n  case,  yet  possesses  certain  unique 
features  themselves  worthy  of  investigation.  Accordingly, 
the  most  extensive  series  of  computations  were  performed 
for  n  =  1. 

On  the  basis  of  both  these  preliminary  calculations 
and  physical  considerations  the  range  of  A  was  chosen 
to  be  [0.000,  0.200]  and  the  range  of  P  was  taken  to  be 
[0.001,  1.000].   Calculations  were  made  for  each  of  the  150 
r-A  pairs;  these  pairs  uniformly  spanned  the  combined  ranges. 
In  a  later  series  of  calculations  r  was  allowed  to  take  values 
as  high  10.00.   In  subsequent  discussion,  the  quantity  c  , 
through  which  A  enters  the  equation  for  k,  will  be  used 
rather  than  A;  however,  r  continues  in  use  since  c  -  -r 
Initial  studies  also  indicated  that  an  evenly  spaced  100-point 
grid  spanning  the  range  [1.0,  11.0]  in  the  radial  coordinate 
was  sufficient  to  study  the  variation  of  the  calculated 
quantities  throughout  the  disk. 

The  dependence  of  the  solutions  to  the  k  equation 
on  cr  r,  and  u  will  be  the  subject  of  the  next  three 
subsections.  One  important  result  should  be  mentioned  at 
the  outset— the  roots  of  the  k  equation  under  almost  all 
circumstances  studied  are  real. 


The  Effect  of  a 

Rather  than  k,  the  absolute  value  of  the  period 


Figure  4   The  quantities  |PER,  |  and  ]PER2| 
are  plotted  agaxnst  radial 
coordinate  for  two  different  values  of 
c, .   Curves  a  and  b  are  |PER2|  and 
| PER,  |  ,  respectively,  when  c,  =  0.361; 
c  arid  d,  when  c,  =  2.709.  In  all 
cases  r  =  0.0757 
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Figure  5   The  dependence  of  IperJ  and 
|  PER-i  1  on  w  is  shown  for  two 
different  values  of  c,.  Curves 
a  and  b  are  | PER2 1 and   | PEI^  | , 
respectively,  when  Ci  =  18.  061; 
c  and  d,  when  c^  =  36.122.   In 
all  cases  r  =  0.075. 
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Figure  6   The  dependence  of  | PER^ |  on 
w  is  shown  for  four  values  c^ 
In  all  cases  r  =  0.075. 
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T 

|  PER  |  =  — - 
2ir|k| 

will  be  the  basis  of  the  following  discussion.  At  each 
point  in  the  disk  there  are  three  k  values,  and  hence  three 
periods.  These  will  be  ordered  and  numbered  according  to 
their  absolute  value—that  is,  IPER.J  <|PER2|  <  |PER3|  .  The 
k  s  will  take  their  numbering  from  the  periods,  for  instance, 
ic1  is  the  value  for  k  which  gives   |  PERX  |  . 

In  Figures  4  and  5,  the  variation  of  IperJ  and 
| PER  (throughout  the  disk  is  plotted  for  four  different 
values  of  c^   The  contribution  of  c±   to  the  coefficients 
of  the  k  equation  are  small  for  small  values  of  c1#  but  they 
are  not  small  for  large  values  of  cr  Yet  the  curves  for 
| PER, |  and | PER2 |  as  a  function  of  radial  coordinate  remain 
essentially  unchanged  from  the  lowest  value  of  c1  to  the 
highest.  Further,  in  each  case,  these  curves  are  not 
sensibly  different  from 

3/2 
| PER |  =  constant  xu 

over  most  of  the  disk  (u  >  2.0)  .   The  constant  is  approximately 
22h  for  IperJ  and  48h  for  |PER2|. 

The  behavior  of  |PER3|  is  different.   In  Figure  6 
|PER3|  is  different.   In  Figure  6  |PER3|  is  plotted  against 
the  radial  coordinate  for  four  values  of  a^.      First,  for 
all  values  of  c,,  the  value  of  IPERj  is  much  greater  than 
the  corresponding  value  of  IperJ  or  |PER2|.  Also  |PER3| 
shows  an  infinity  (k"3  goes  to  zero)  whose  location  is 


PER, I  (curves 


(curves  a  and 


Figure  7   The  dependence  of 
b  and  d)  and  | PER2 
c)  on  a)  for  two  values  of  r. 
Curves  a  and  b  are  for  the  case 
T  =  0.001;  curves  c  and  d  for  the 
case  r  =  0.01.   In  both  cases 
C;L  =  0.361. 
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Figure  8   The  dependence  of  IPERjJ  (curves 
b  and  d)  and   |PER2|  (curves  a 

and  c)  on  a)  for  two  values  of  *■■ 

T  .  Curves  a  and  b  are  for  the 
case  T  -  0.10,  and  curves  c  and 
d  for  the  case  r  =  1.00.   In  both 


cases  c1  =  0.361. 
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2.0  4.0  6.0  8.0  10.0 

Radius  (Stellar  Radii) 


Figure  9  The  dependence  of  |PER3|  on 
w  is  shown  for  two  values  of 
T  :  r  =0.001  and  r  =  0.01. 
In  both  cases  c,  =  0.361. 


92 


2.0 


4.0  6.0  8.0 

Radius  (Stellar  Radii) 


10.0 


93 


dependent  on  c...     Beyond  this  infinity,  further  out  in 

the  disk,  |PER3|  falls,  off  first  rapidly  and  then  slowly 

5 
with  increasing  w.   It  approaches  a  value  of  about  1.8  x  10  h 

near  the  outer  region  of  the  disk. 


The  Effect  of  A 

There  are  many  similarities  with  the  preceding 
case.   The  dependence  of  IPER.J  and  |PER2I  on  the  radial 
coordinate  for  four  different  r  values  is  shown  in  Figures 
7  and  8.   The  curves  change  only  a  little  with  large  changes 
in  r.  Only  for  r  =  1.00  (curve  c  of  Figure  8)  do  any  of  the 

curves — and  then  only  |PER2I —  show  a  noticeable  departure 

3/2 
from  u  '      behavxor. 

The  quantity   |PER,|is  much  more  dependent  on  r. 

Examination  of  Figures  9  and  10  shows  that  as  the  value  of 

r  increases  the  singularity  in  |PER3| moves  inward;  by 

r  =  0.25  (curve  c  in  Figure  10),  it  is  interior  to  the  first 

grid  point.   Further  there  is  a  general  decrease  of  |PER3| — 

away  from  the  singularity — as  r  increases.   For  r  =  1.00 

(curve  e  in  Figure  10),  |PER,|  reaches  its  maximum* 300 

days,  at  u  =  7.6  and  decreases  outward. 


Extreme  Cases 

In  Figure  11,  the  radial  dependence  of  the  three 
periods  when  both  c1  and  r  are  large  is  displayed.   In  the 
outer  part  of  the  disk,  |PER2|  shows  the  same  departure  from 


Figure  10   The  dependence  of  |PER  |  on  u  is 
shown  for  three  values  of  T: 
r  =  0.050  (c),  r  =  0.250  (d) ,  and 
T  =  1.00  (e) .  Again,  c,  =  0.361 
in  all  cases. 
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Figure  11   The  dependence  of  the  three  periods 
on  w  is  shown  when  c-i  =  36.122  and 
r  =1.00.   |PERX|  is  a,  |PER9|  is  b, 
and|PER3|  is  c.  z 
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Figure  12   The  dependence  of  |PER,|  (a), 
|  PER,  |  (b)  ,  and  |PER3r(c)  on 
id  when  c,  =  0.00  and  r  =  0.00 
is  shown. 
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3/2 
w  '   behavior  that  is  found  for   | PER2 |  in  curve  c  in 

Figure  8.  Otherwise  [PERjJ  and  |PER2|  show  the  same  insensi- 
tivity  to  the  exact  value  of  c,  or  r that  was  shown  in 
earlier  cases.   The  location  of  the  singularity  in  |PER  | 
is  the  result  of  competing  trends  for  large  c..  and  large  r. 
Large  values  of  c1  move  the  singularity  outward,  while 
large  values  of  T  move  it  inward  toward  the  star.   The 
location  of  the  singularity  is  approximately  midway  between 
its  location  for  large  c1   alone  (o,  =  36.122  in  Figure  6) 
and  the  position  of  the  singularity  when  only  r  is  large 
(T  -  1.00,  curve  e  in  Figure  10).   In  the  outer  regions 
of  the  disk,  one  sees  effects  very  much  like  those  in 
Figure  10. 

The  radial  dependence  of  these  quantities  for  small 


PER 


2 


parameter  values  is  shown  in  Figure  12. JPER1  |  and 
(a  and  b  in  the  figure)  confirm  the  pattern  of  relative 
insensitivity  of  the  two  quantities  to  changes  in  the 
parameters.   |PER3|  (curve  c)  is  quite  large,  shows  no 
singularity,  and  increases  in  the  disk  as  co1/2. 

The  parameter  r  determines  the  presence  or  absence 
of  imaginary  roots.   Calculations  for  large  values  show 
that  imaginary  roots  appear  in  the>  outer  regions  of  the 
disk  at  about  r  =  2.50.  As  r  increases  beyond  this  value 
larger  and  larger  regions  of  the  disk  exhibit  imaginary 
roots.   The  onset  of  imaginary  roots  is  insensitive  to  c.  . 


Figure  13   The  dependence  of  the  three  periods 

on  a)  is  shown  for  three  values  of  n — 3 
(a) ,  26  (b) ,  and  41  (c) .  In  each  case 
c1   =  0.361  and  T  -  0.075. 
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Figure  14   The  quantity   |PER,|  at  oi  ■  .1.1 
is  plotted  against  n.  For  all  n 
c^  =  0.361  and  r  =  0.0750.   The 
solid  curve  is  52.0/n. 
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Figure  15   These  curves  illustrate  the  effect 
of  extreme  values  of  c.  and  r 
for  n  =  41.  Curve  a  corresponds 
to   I  PER, |  when  c±   =  0.000  and 
r  =  0.001.  The  other  two  periods 
which  fall  quite  close  to  this 
curve  are  not  shown.  The  curves 
b  (| per 3|),  c  (|PER  |),  and  d 
(|PER \y   are  for  thS  case  in  which 
c1  =  36.122  and  r  =  1.00. 
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The  Case  n  >  1 


General  Features 


The  appearance  of  a  singularity  in  |PER3|  is 
unique  to  n  =  1.   Figure  13  shows  all  three  periods  plotted 
as  a  function  of  the  radial  coordinate  for  three  values  of 
n — 4,  26,  and  41  (curves  a,  b,  and  c  respectively).  Again 
one  sees  the  a>3//2  dependence  in  |  PEI^  |  and  |  PER2 1  ,  and  now 
in  | PER.,  |  as  well.  Note  that  the  separation  among  the  three 
periods  decreases  with  increasing  n.  More  importantly,  the 
absolute  values  of  the  periods  fall  off  markedly  from  n  =  4 
to  n  =  26,  and  less  markedly  from  n  =  26  to  n  ■  41. 

This  behavior,  a  decrease  with  increasing  n,  is 
confirmed  in  Figure  14.   Here  the  value  of  |PER.jJ  at 
u  =  1.1  has  been  plotted  against  n.   The  solid  line  in  the 
figure  is 

Iper  I       -  0-520  x  10+2  . 
lPERl'  w=l.l  ~ 

the  fit  for  n  greater  than  10  is  remarkable.   The  other 
two  periods  show  similar  behavior. 

Large  n  and  the  Effect  of  c.  and  r 

Figure  15  shows  the  effect  of  extreme  values  of  the 
constant  c,  and  r  on  the  three  periods  in  the  case  n  =  41. 
In  general,  these  and  other  calculation  show  some  dependence 
on  r,  and  almost  none  on  c...  Further,  for  those  n  values 
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tested,  no  roots  of  the  k  equation  are  imaginary,  even  for 
r  >  5.00.   It  would  be  prohibitively  expensive  to  perform 
the  extensive  calculations  done  for  the  n  =  1  case  on  very 
many  n  values.  But  one  can  find  at  least  qualitative 
support  for  these  findings  in  the  k~  equation. 
Writing 

k  =  n  A 

and  substituting  this  into  the  equation,  one  finds  that 

c  2  c 

A3   +   3f   A*  +    (  -|  +  3f2  -  ^j  -  -5^  )    A 

a)  nun 


c9f  ,        2c   f  c  f  f3 

+      {      2  22        ,   52  2J 

u  n  a)  2uJn         n 


-2 
The  integer  n  appears  in  five  terms  as  n   .   For  large  n  these 

terms  are  negligible  and  c,  occurs  only  in  these  terms.   For 
large  n,  at  least,  the  solutions  of  the  k  equation  should  not 
depend  on  c,. 

For  large  n  the  equation  for  A  is 

A3  +  3fA2.  +(-  4  +  3fV  +  (-1!  +  f3')  =  °   • 
For  r  in  the  range  [0.001,  1.00]  and  w  in  tl.0,  ,11.0], 

|-*-|   <   0.267   . 

f   (A) 

The  terms  which  contain  r  are  never  dominant;  generally  they 
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may,  to  a  first  approximation,  be  neglected.  If  they  are, 
one  has  that 


A3  +  3fA2  +  3f2A  +  f3  -      0  . 


Clearly, 

A  a  -f   . 
This  gives 

2ttT    T  u3/2 
IPERl  =  °  =  _2 


nf      n 


=  48^5  w3/2 


in  hours.   At  a)  =1.1 


|PER|   ,  ,  -  — 
1    ' w=l . 1   n 


This  result  compares  favorably  with  the  curve  in  Figure  14. 

The  approximate  nature  of  this  result  should  be 
emphasized.  Calculations  show  that  in  all  cases  there  are 
three  different  roots.  When  r  is  large,  one  of  these  roots 

may  be  expected  to  depart  from  an  id'   radial  dependence, 

2 

becoming  more  like  fi>,  . 
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Sample  Flow  Variable  Solutions 

As  was  discussed  in  Chapter  IV,  there  can  be.  no  fully 
determined  solutions  from  this  analysis.   These  solutions 
can  be  determined  only  within  a  factor  S  r"(w)  ,  which  must  be 
given.   The  qualitative  nature  of  the  solutions  can  be 
studied  by  choosing  some  simple  S  r-   and  calculating  the 
solutions  as  functions  of  u),  <J>,  and  t.  A  simple  choice  is 


Snk(w)  "  IT"  6nl  6kk\ 
lOw         1 


where   6  ,  and  Sj-j-  are  Kronecker  deltas.   Figure  16  shows 
u,^   and  u,,  as  functions  of w near  the  star  [1.00,  3.00]  and 
for  three  different  times:  0  hours,  40  hours,  and  80  hours. 
Figure  17  shows  the  quantity  [u,,  (90°)  -  u,,  (270°)] 
at  two  times:  40  hours  and  80  hours.   On  the  basis  of  this 
figure  one  would  expect  a  variation  in  the  emission-peak 
separation  (See  Figure  1} . 

Several  general  features  stand  out  in  this  case 
which  were  also  present  in  other  sample  calculations  not 
shown  here.  First,  all  the  cases  are  complex.  Observations 
of  phenomena  which  occur  in  the  disk  and  which  reflect 
these  solutions  will  show  extremely  complex  temporal  behavior, 
Second,  the  density  for  the  above  S  j-  corresponds  to  a  change 
in  density  which  is  never  greater  than  1/10  of  the  steady- 
state  value.   Yet  the  velocities  approach  1/3  of  the  charac- 
teristic steady-state  velocity — near  the  practical  limit  at 
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which  the  linearization  procedure  can  be  considered  valid. 
This  suggests  that  only  small  changes  in  density  are 
associated  with  appreciable  velocity  changes.  Finally,  near 
the  stellar  surface  the  velocities  of  neighboring  portions 
of  the  disk  may  be  quite  different  and  the  relative  velocity 
between  two  such  points  may  change  rapidly. 


Comparison  with  Observation 
Actual  Stellar  Values  for  Vg 

The  quantity  A  always  appears  in  the  equation  for  k 
in  the  form 

dl — 2   * 

2FN     V 

s 

2 
The  quantity  Vg  ,  the  square  of  the  equatorial  rotational 

velocity  of  the  central  star,  appears  in  the  definition 
of  r.   Variations  in  these  quantities  have  been  viewed  as 
changes  in  the  physical  conditions.  As  far  as  the  mathematics 
is  concerned,  however,  these  variations  could  be  the  result 
of  changing  V  . 

Actual  stellar  values  for  V  lie  in  the  range  from 
the  present  value  (half  the  break-up  velocity)  to  the  break- 
up velocity.   This  range  can  give  at  most  a  factor  of  4 
decrease  in  c±   and  r.   This  value  is  small  in  comparison 
with  the  range  of  values  of  c1  and  A  over  which  calculations 
were  performed.   Considering  the  results  of  this  chapter. 
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it  is  clear  that  the  calculations  are  very  little  influenced 
by  the  exact  value  of  V  .   It  is  also  clear  that  calculations 
do  not  have  to  be  performed  for  a  variety  of  Vg  values.  Vg 
does  appear  in  the  dimensioning  parameters.   For  a  value 

of  V  other  than  250  km/s,  the  periods  should  be  reduced  by 

s 

the  factor  V  /250,  where  V  is  in  km/s. 
s  s 


Qualitative  Aspects  of  the  Predicted  Spectral  Variations 

Even  in  the  absence  of  fully  determined  solutions, 
one  can  determine  many  of  the  qualitative  aspects  of  the 
spectral  variations  which  these  solutions  predict.  A 
knowledge  of  these  qualitative  aspects  is  adequate  to 
discuss  much  of  the  observational  material. 

In  Chapter  I,  the  spectral  variations  associated 
with  Be  stars  were  divided  into  three  rough  classes  according 
to  the  characteristic  time  of  the  variation;  these  were: 

1)  long  term  (>  10  years) , 

2)  medium  term  (1  to  10  years) , 

3)  short  term  (<   1  year) . 

Clearly  the  phenomena  which  may  result  from  these  solutions 
fall  into  the  last  category,  short  term  variations  (the  only 
exception  being  the  case  for  n  =  1  and  k  =  k~3)  .  Depending 
on  the  value  of  n,  the  periods  in  the  inner  part  of  the  disk 
range  from  several  days  to  fractions  of  a  day,  and  for  larger 
n  values ,  fractions  of  an  hour. 

For  a  given  n,  the  values  of  the  periods  change 
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throughout  the  disk,  and  the  presence  of  several  n  values 

at  once  cannot  be  precluded.   The  multiplicity  of  periods 

suggests  that  phenomena  dependent  on  these  solutions  should 

show  a  complex  variation  with  time.  Although  one  period 

may  dominate  the  observations ,  in  general  no  simple  periodicities 

should  be  expected. 

Following  the  taxonomy  of  McLaughlin  (1961) ,  the 
two  major  types  of  spectral  change  which  could  be  expected 
to  result  from  subscript-one  solutions  are  V/R  variations 
and  E/C  variations.   The  total  emission  in  a  line  in  the 
circumstellar  envelope  should  be  more  dependent  on  the 
density  than  on  the  fluid  velocity.   The  ratio  of  the 
intensity  of  the  violet  emission  lobe  to  that  of  the  red 
lobe  is  strongly  dependent  on  the  fluid  velocity — particularly 
the  radial  component.   In  the  sample  solutions  which  were 
calculated,  the  amplitude  of  the  density  variation  is  much 
smaller  than  the  amplitude  of  the  velocity  variation.   The 
presence  of  these  solutions  should  be-most  evident  in  the 
V/R  variations. 

The  largest  period  for  the  case  n  =  1  is  an  anomalous 
one.  For  typical  parameter  values  this  period  is  greater 
than  10  years  over  most  of  the  disk.   Long  term  changes 
observed  in  Be  stars  have  generally  been  large  changes — not 
ones  which  one  would  expect  for  solutions  to  the  linearized 
equations .  For  certain  combinations  of  the  parameters , 
| PER- |  for  n  =  1  is  for  most  of  the  disk  on  the  order  of 
years  and — for  large  enough  r — fractions  of  a  year.  Further, 
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for  large  values  of  r  |PER3|  changes  only  slowly  with  u 
(see  curve  e  of  Figure  10) .   This  suggests  that  any  observable 
effect  of  |PER3|  for  large  r  will  possess  a  simple  periodicity. 

It  is  interesting  to  compare  this  with  the  81-day 
period  found  in  HD41335  by  Mrs.  Peters  (1972).   The  period 
is  based  on  the  appearance  and  absence  of  emission  line 
structures  in  spectra  taken  at  a  variety  of  times.   If  one 
uses  suggested  stellar  rotational  velocity  of  585  km/s, 
the  maximum  in  curve  e  of  Figure  10  corresponds  to  a  period 

of  slightly  over  100  days.   With  somewhat  larger  values  of  r 

the  agreement  would  be  better. 

A  physical  reason  for  the  large  r  value  must  be 

provided  before  a  strong  case  can  be  made  for  this  explanation. 

Mrs.  Peters  felt  the  periodicity  was  the  result  of  a  companion; 

Periodicity  phenomena  in  other  cases  have  been  explained  as 

effects  of  a  companion,  but  the  binary  nature  of  HD41335  has 

not  yet  been  demonstrated. 


Comparison  with  Observations  Reported  in  the  Literature 

Variations  in  spectral  features  of  Be  stars  which 
become  evident  in  a  year  or  less  have  been  known  to  exist 
for  some  time, yet  no  large  body  of  observations  on  them 
exists.   This  is  particularly  true  for  variations  over  days, 
fractions  of  a  day,  and  fractions  of  an  hour — this  last 
requires  "state  of  the  art"  observing  techniques  (Hutchings 
et_al. ,  1971).   The  principal  studies  in  this  time  region 
(several  days  and  shorter)  will  be  discussed  below.   There 
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are  two  reports  in  the  literature  of  strictly  periodic 
variations;  both  are  V/R  variations.   Lacoarret  (1965) 
found  a  seven-day  periodicity  in  HD174237.   Hutchings  (1970) 
found  a  0.7-day  periodicity  in  Y  Cas.   In  the  second  case 
(  y  Cas) ,  where  there  was  a  20-minute  time  resolution, 
numerous  secondary  fluctuations  partially  obscured  the  0.7-day 
periodicity.   There  was  also  a  variation  in  the  peak  separation 
of  the  double  emission  profiles  with  the  same  period.   The 
first  case  could  be  explained  by  solutions  with  small  n 
values  and  with  appreciable  amplitude  only  in  the  outer 
parts  of  the  disk.  One  would,  however,  expect  some  evidence 
of  other  periodicities;  these  are  not  present.  The  second 
case,  with  a  shorter  period  and  with  many  apparently  irregular 
fluctuations  clearly  present,  is  qualitatively  a  kind  of 
behavior  these  solutions  might  predict. 

The  observations  most  relevant  to  the  work  of  this 
chapter  are  the  additional  observations  of  y  Cas  by 
Slettebak  (1969)  and  Hutchings  (1969)  and  of  HD142926  as 
well  as  k   Dra  (Hutchings  et  al. ,  1971) .   These  observations 
are  characterized  by  the  following: 

1)  highly  complex  time  dependence, 

2)  rapid  temporal  variations, 

3)  V/R  variations  in  a  complex  emission-absorption 
structure . 

In  the  observations  of  HD142926  and  *  Dra  temporal  resolution 
was  down  to  2  minutes;  variations  were  observed  to  this  limit. 
Times  of  variation  are  not  periods,  but  the  time  over  which 
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a  change  was  noted.   These  compare,  not  to  periods,  but 
rather  to  a  number  like  |per|/4  or  |per|/6.  If  large  n 
values  (approximately  to  100)  are  included,  the  predictions 
of  the  analysis  and  the  observations  agree  on  all  three 
characteristics . 


CHAPTER  VI 
EQUATION  FOR  ENERGY  TRANSPORT 

Form  of  the  Equation  of  Energy  Transport 

The  Distinctive  Nature  of  the  Energy  Equation 

In  one  respect,  the  equation  for  energy  transport 
stands  in  sharp  distinction  to  the  other  two  hydronynamical 
equations;  the  very  formulation  of  a  general  equation  for 
energy  transport  requires  an  intimate  knowledge  of  the 
interaction  between  the  radiation  field  and  the  disk  fluid. 
Absorption  and  emission  in  the  disk  determine  the  extent 
of   ionization  throughout  the  disk,  and  since  these  processes 
act  as  heating  and  cooling  mechanisms,  the  temperature.  Yet 
radiation  terms  appear  not  at  all  in  the  continuity 
equation,  and  as  small  effects  (See  Chapter  II)  in  the 
equation  of  motion. 

No  effective  treatment  of  the  radiation-fluid 
interaction  which  is  applicable  to  the  envelopes  of  Be  stars 
exists.  At  the  onset,  then,  the  form  of  the  equation  for 
energy  transport  must  be  more  approximate  than  the  initial 
formulation  of  either  of  the  other  two  flow  equations. 
For  this  reason  and  because  of  the  mathematical  complexity 
which  marks  even  the  simplest  form  of  the  energy  equation, 
the  equation  for  energy  transport  merits  a  separate  discussion, 
In  the  absence  of  an  effective  treatment  of  the 
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radiation-fluid  interaction,  a  minimum  condition  can  be 
imposed.  One  can  demand  formal  similarity  to  the  analysis 
of  Chapter  III  and  IV.  A  simple  criterion  provides  this. 
Given  a  candidate  equation,  the  following  procedures  are 

followed: 

(1)  typical  disk  conditions  are  imposed; 

(2)  each  flow  parameter  is  expressed  as  the  sum  of 
two  terms — a  steady-state  axially-symmetric 
part  (subscript- zero  component)  and  a  part 
containing  the  angular  and  temporal  dependence 
(subscript-one  component) ; 

(3)  the  equation  is  linearized,  and  a  relation 
which  must  be  satisfied  by  the  subscript-zero 
terms  is  constructed. 

The  criterion  is  that  the  static  isothermal  solutions  to 
the  equations  of  continuity  and  motion  also  be  solutions 
of  the  relation  found  in  step  (3)  above.  The  equation 
for  the  subscript-one  terms  may  then  be  examined  in  the 
context  of  previous  work. 

In  the  remainder  of  this  chapter,  an  equation  for 
energy  transport,  the  most  general  one  possible  in  the 
absence  of  an  effective  treatment  of  the  radiation-fluid 
interaction,  will  be  proposed.   The  procedures  and  criterion 
discussed  in  the  previous  paragraph  will  be  followed.   The 
mathematical  aspects  of  the  analysis  followed  in  Chapter  V 
can  be  applied  then  to  the  linearized  equation  for  the 
subscript-one  terms  originating  in  the  energy  equation. 
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The  final  sections  of  this  chapter  examine  first  the 
agreement  with  the  work  of  Chapters  III,  IV,  and  V  and  then 
the  adequacy  of  the  approximate  form  of  the  energy  equation. 


An  Approximate  Form 

Any  radiation  effects  enter  only  indirectly  through 
the  internal  energy  and  possibly  through  a  somewhat  stellar 
mass  term  (originating  in  the  equation  of  motion  and  included 
here  only  for  completeness) .  A  simple  three-element  fluid 
consisting  of  protons,  electrons,  and  neutral  helium  will 
be  used.  The  equation  for  energy  transport  could  proceed 
in  the  same  manner  as  the  equation  of  motion;  that  is, 
the  equation  for  energy  transport  could  be  derived  from  an 
integral  involving  a  conserved  quantity — the  thermal  energy. 
The  equation  given  below  is  an  equation  for  energy  transport 
in  the  presence  of  viscosity.   It  is  from  the  mathematical 
standpoint  the  same  order  of  approximation  as  the  Navier- 
Stokes  equation.   In  practice,  the  equation  is  derived 
heuristically  (Landau  and  Lifshitz,  1960) .   The  equation  is 

2  2 


9   C-fr  +  ^c    +  *cpc  )  «  -v  •  [pcv-c  (  -2  +  wc  +  #J 

c 


at    2    °  °        c  c  ;      u  c  c   2   c   c 


-  ^c  '  ac  "  KcVTc  J 
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In  this  equation  the  symbols  not  used  before  are: 

(1)  e  for  internal  energy  density,  . 

(2)  w  for  enthalpy, 

(3)  K     for  thermal  conductivity, 

c 

(4)  a1   for  the  viscosity  stress  tensor, 

c 

(5)  $  for  the  gravitational  potential. 

c 

The  product  v"  •  a  *  has  elements 


(v.  •a').=v.a' 
vvi   uc7d    x  c,x,d 


where 


3v0     3v    ,      3v„      3V, 
C'1'^      Sx,     3Xjl   3   ij  3x£      3x£   1:J 


D 


and  repeated  indices  indicate  a  sum.  Note  that  the  subscript 
c  which  would  normally  appear  with  v  and  x  above  has  been 
supressed  in  the  expression  to  avoid  confusion  with  the 
indices  of  summation.   The  internal  energy  density  and 
enthalpy  for  the  three- element  gas  are 


3kTc  .  XI 

e  =  +  — 

c   2ym    pm 


and 

5kTc  -  XI 

W   =  — : +  — - 

c   2ym  um 
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respectively,  where  X  is  the  mass  fraction  of  hydrogen 
and  I  is  the  ionization  potential  of  hydrogen. 


Relative  importance  of  the  Terms  of  the  Equations 

A  new  quantity,  the  thermal  conductivity,  appears 
in  the  energy  equation.   Since  the  mass  of  the  electron  is 
so  much  smaller  than  that  of  either  the  proton  or  the  helium 
atom,  the  thermal  conductivity  in  the  three-species  mixture 
is  essentially  determined  by  the  electron.   The  thermal 
conductivity  for  the  electron  component  of  the  gas  is 


Kc  = 


1.952  x  10~4  T^/2  3\ 
c ±_ 

log  Ud)  Fx 


where  6*.  and  6,  are  dimensionless  numbers  of  order  unity 
given  by  Spitzer  (1962) .   The  ideal  gas  law  will  be  used  to 
express  the  last  term  •  in  the  equation  as  a  function  of 
pressure  and  density. 

The  dimensionless  variables  discussed  in  Chapter 
III  are  introduced  into  the  equation  for  energy  transport 
giving 


3t   z  *  ymVZ     V*R 

o      °  ° 
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->•  •*       M  G  ±  "  

,  5PV  .  XI PV  _  (        s  V£  _  __TJ /   a.  } 

+  ~T~  ~ — 2    «•  H2      r   II  V  R    i  oj'i 

2  unfiT      V~R     r   Jlovo*b     UJ   J 


<vl      '   7/2  ,  „  ,  5/2 


'  o 


7/2   _   5/2     D   "I 


where 


1 


cr',  ,  =  =•  a' 


and 

1.952  x  lO"4^ 

K    =    ^^ 

Both  sides  of  this  equation  have  been  divided  by  the  quantity 

(II  V2/T  ).   The  quantity  a.,   is  of  order  unity. 
x  o  o  o  ID 

The  coefficients  should  be  measures  of  the  relative 
importance  of  their  respective  terms.   To  evaluate  them  the 
characteristic  values  (following  Chapter  III)  are  taken 
to  be: 


(i)  no  =  pd  , 


(2)  RQ  =  Rg   , 

(3)  V  =  V  /2tt 

o    ** 
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<4>   To  =  VVo   ' 
<5>   Po=\noVo   • 


Table  1  contains  the  values  used  for  the  quantities  Rg,  Vg/ 
and  p, .  The  equation  is  rewritten 


ft  (  Al  ^     +   A2  P  +  A3P  +  A4  F  )  =  "V  "   [BlpV^ 


+  B2PV 


+  B3  pv  +  B4 


#+B5,v.?ij)+B6  (|)5/2V(|)] 


The  values  for  these  coefficients  are  given  in  Table  2. 
Two  terms  merit  special  comment.   These  are  the 
last  two  terms,  which  are  much  smaller  than  the  other 
elements  of  the  equation — so  small  that  they  may  be  dropped. 
In  the  case  of  the  last  term,  (Bg) ,  which  is  clearly  non- 
linear, this  is  fortunate.  Even  though  both  are  small,  the 
disparity  between  the  two  values  B5  and  Bg  is  at  first 
thought  disturbing.   The  coefficient  B5  is  the  inverse  of 
the  Reynolds  number;  the  thermal  conductivity  term  is 
approximate ly 


KVI    «  . V2     k„ 


k  J        '  «  TT  ,7    l  k 


151(4)      -r=*r(f)« 


IIR     P   -  *  RIIV    *     n  R  V  C      PN 

o  o  o  o  o        o  o  o  p 

where  C  is  the  specific  heat  at  constant  pressure  and  PN  is 
P 
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TABLE  2 

The  Values  of  the  Coefficients  Appearing 
in  the  Energy  Equation 


Value 
1.00 

1.50 

0.88 


Coefficient 

Al 

A2 

A3 

A4  . . 

Bl 

B2 

B3 

B4 

B5 

-1.22  x  10+2 


0.50 
2.50 
0.88 


-1.22  x  10+2 


8.11  x  10"13  (protons) 
or  8.81  x  10-11  (helium) 


Bc  2.02  x  10"6 

o 
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the  Peclet  number.  Now  the  ratio 


5S.  »  i 
pn  >y    1  ' 


yet  the  conventional  wisdom  of  hydrodynamics  is  that  these 
two  numbers  should  be  roughly  equal.   This  discrepancy 
appears  to  be,  not  an  error  in  the  analysis,  but  rather  a 
peculiarity  of  the  ternary  mixture  chosen  for  study  here. 
The  principal  agent  for  thermal  conduction  is  the  electron 
component  of  the  gas  (Chapman  and  Cowling,  1970) ,  while 
protons  and  neutral  helium  provide  most  of  the  viscosity. 


A  Linearized  Equation 
Separation  of  the  Time  and  Angle  Dependence 

The  equation  for  energy  conservation  is 

il  [-f*1  ♦  5  ♦ »'  --fr©]  --'  -[p*(4  +  §  + 


The  flow  variables  are  separated  into  subscript-one  and 
subscript-zero  quantities  as  before.   These  are  then 
substituted  into  the  energy  equation  which  becomes 

it  ft  <*o  •  5o  +  2*o  •  ?1>  +  ^  <V  *o>  +  Bpo  +  Bpl 
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+  ~2~       ^^       FN      L  r   J     .  FN      >     r  J  J  .  2      uvo     V 


♦  2^.^)1+^  <?0 -  v0>  *-;2sa  <*o  •  ?o>  +  Bpo?o 


+  B 


-    5  _  i    5PoV1   5plvo 


povl  +  BplVo  +  2  Povo  +  —2—  + 


2 


1   #± 

FNr 


(vo  po  +  vo  pl  +  vl  po)  J  • 


The  terms  which  involve  subscript-zero  quantities 
alone  are  collected  on  each  side  of  the  equation  and  are 
set  equal  to  each  other.   The  equation  which  must  be  satisfied 
by  the  subscript  zero  solutions  is 


^[%<vo-vo.  +Bpo  +  ^-^|  ]  - 


V  '[  povo(vo  •  vo)  +  Bpovo  +  TT  vo  -  FN?  povoJ 


Since  the  subscript-zero  quantities  are  first  steady-state 
and  axially-symmetric  (there  is  not  time  or  angle  dependence) , 
the  equation  for  the  subscript- zero  terms  becomes 


1   3   f  "poVa)Q  v2  .   _  w    .  ^o^o  .  MpoVcoo  1 
Jo  tt     L  -1 VD  +  wBpoVa)0  +  2 FNr-  -I 
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|..[^  ♦-?*'-♦' 3^-'*  P0vzo]  -•« 


where 


2  =  v2  +  v2  +  v2 
o    wo    zo    q>o 


The  subscript-zero  solutions  used  here  are  also 
static  and  isothermal  solutions  in  which 


=  v(u)  *# 


These  solutions  satisfy  the  equation  above  identically.   This 
form  of  the  energy  equation  satisfies  the  criterion  proposed 
in  the  first  section  of  this  chapter. 

An  Equation  for  the  Subscript-One  Terms 

The  static  isothermal  solutions  do  satisfy  the 
relation  for  subscript-zero  quantities  which  results  from 
the  energy  equation.  Accordingly,  the  static  isothermal 
solutions  used  for  subscript-zero  quantities  in  Chapters 
IV  and  V  will  be  used  here.   They  are: 

<«   vu)o  =  °  ' 
2tt 

<2>  %o  ■  -TT2    ' 
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(3)  vzo=0. 

(4)  P0-«p[-|(i-'i)]. 

(5)  p0  =  *exp[-i  (1    -     |)]  • 
The  equation  for  the  subscript-one  terms  is 

j|  t  Pov*ov*l  +  I  v*o  pl  +  Bpl  +  41  "  fs£  '•  = 


2 


-(i»[^MP0^^    PJ     vul  +  t+[(^  + 


3 

B(;o  +  T2  "  Sfif  '  %i  +   (  T2  "  ?HF  V  +  Bvtoi   "l  +  !  v»opl  J 


„      r        P    iovf  5(i)p  um  -i 


■;i}  [(f^b*o*»'o  +  Ta-l»l)  V 


♦  c  3*  -  ite  +  bv*o  )  pi  + 1  wi  ] 

8?     L   I       2  po         2  FNr  J      zl  J 
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The  subscript-zero  quantities  are  known  functions, 
and  consequently  all  derivatives  of  these  known  quantities 
may  be  obtained.   Further,  the  mathematical  form  of  these 
functions  and  their  derivatives  suggests  that  the  substitution 
of  expressions  valid  for  z/w  <  1  will  simplify  the  calcula- 
tions at  no  great  cost  in  usefulness.   The  new  dependent 
variable 


vi  ■  vi  exp  c  -=4- ) 


is  introduced  into  the  equation.  The  result  of  all  these 
steps  is 

3vd>l     m'2*2   3Pi    ■*    3Pi      3p,     ,   3p, 

3t       2   at   2  at  a  "TE  fnF  It  " 

-f  B  +  OT  1  ve  _f  if|2       jf   ^      _     ve 

l  a)  +  2  J  vwi  I  —2~  +  B  +  —  J  -air  +  u  a^  I  "fnF  i 

f  3o)f2  .  B  ■  51   _    1   .  3v»l 
1   2    w   2w   FNru  J   3(f> 

v   2         FNr  }      3<J>    2   3<J)  +  "FN   3z  I  ~r"  J 
1  2  +  B  +  T  J  -17"  ' 


In  the  rest  of  this  chapter  the  superscript  e  will  be  dropped; 
the  velocity  will  contain  an  absorbed  exponential  unless  other- 
wise noted. 
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The  mathematical  steps  carried  out  thus  far  closely 
parallel  the  analysis  of  the  equations  of  continuity  and 
motion  performed  in  Chapter  VI.   The  set  of  four  equations 
in  that  chapter  which  correspond  to  the  same  stage  of 
analysis  as  the  equation  above  are  simpler  mathematically 
than  this  one  equation.   This  mathematical  complexity  can 
be  an  important  limitation  to  the  usefulness  of  the  energy 
equation. 

Two  questions  suggest  themselves.   In  Chapter  IV  a 
series  of  approximations  tailored  to  the  physical  problem 
and  to  the  mathematical  peculiarities  of  the  equations  were 
applied  to  the  four  linearized  scalar  equations  representing 
mass  and  momentum  conservation.   These  approximations  resulted 
in  a  set  of  three  solvable  equations  for  three  integrated 
flow  variables.  Are  the  results  of  that  analysis  consistent 
with  the  result  of  the  application  of  the  same  series  of 
approximations  to  the  energy  equation?   If  they  are,  does 
the  equation  for  energy  transport  above  allow  the  construction 
of  less  approximated  solutions?  The  first  question  will  be 
examined  next. 

Relationship  to  Previous  Work 

Reduction  to  Comparable  Form 

At  this  point  the  flow  will  be  assumed  to  possess 

symmetry  about  the  equatorial  plane.   Further,  the  equation 

for  energy  transport  will  be  integrated  with  respect  to  z 

over  the  interval  [-z  .  +  z  ]  where  z  is  a  function  of  u. 

o     o        o 
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Doing  this  will  require  the  following  integrations: 


(1)     I     °    FTdz     ' 
-zo 


"zo 


T  +2      ^(J)! 

/•a\  •       f       O 

(3)    i   L  -w 


°  -*1      dz.     , 


"zo 


+Zo        1  3P1 


o 


<=>  £°  £  ( ¥ ) az 


The  rather  tedious  mathematics  required  to  construct  suitable 
expressions  for  the  integrals  is  contained  in  Appendix  E;  the 
resulting  expressions  for  these  integrals  are: 


«>    ti--£) 


3U 


ail     A' 


.       +   --.'■  ■  u  , 

a)   7^  ;   3(0     ...4    ail 


2u>  w 


m   f  1  -  A'  %  ^1 
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<4>-  ■(*-  -M'ft£   ' 


u   20)3  '■•  **"" 


(5)     Zl  ° 


CO 


The  new  parameters,  A' and  A',  are  defined  by  the  expressions 


/+Zo  ■*%!«■  =  A'/+!°  v  ,    dz 


-z. 


-z      vtol 


and 


C°  2%idz  =  A''-*°  V  dz  - 


respectively,  Both  of  these  parameters  are  considered 
constant . 

One  departure  from  the  pattern  of  analysis  fround 
in  Chapter  IV  is  made  here.  The  quantity 


+z 

Pi  -  /"-■  Pl  dZ 
o 


is  defined.  This  removes  the  need  for  an  assumed  relation 
between  p,  and  p.;  it  is  this  relation  which  defines  the 
parameter  r.  Otherwise  the  steps  outlined  in  the  last 
section  of  Chapter  IV  are  followed  here.  The  result  of  all 
this  is 
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-( I +  IS)   uu  "  ^7  ua  '  (  — +  - +  5s     m„2 


A'      ,  9"*1  ■'        f  o>2f3  +  Bf      .        M_  1   !!l  .  5f  *J 


The  Question  of  Consistency 

The  approximated  form  of  the  energy  equation  and 
the  results  of  Chapter  IV  and  V  are  consistent  if  the  solutions 
worked  out  in  these  chapters  satisfy  that  equation.  To 
test  this,  the  relation 

dPl  =  reff  dal 

is  substituted  into  the  equation  derived  at  the  end  of  the 
last  section.  The  expression  for  the  flow  variables  (u^, 
u   ,  and  a,)  derived  in  Chapter  IV  are  used  to  evaluate 

0)1         1 

terms  in  the  energy  equation  which  contain  these  variables 
or  their  derivatives.  One  then  has  an  algebraic  expression 
which  can  be  used  (See  Appendix  F)  to  develop  an  expression 
for  r  ff.  This  expression  has  the  form 


?   ££  =  £0c(r),n,A#A"JA») 
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The  notation  k(T)  signified  that  r  enters  only  through  k. 
This  equation  depends  directly  on  n-^A^lT1  r  and  A*, 
and  indirectly  on  r.  However,  the  quantities  A*  and  A  have 
little  influence  on  the  equation.  A  simple  procedure 
suggests  itself.   Values  are  chosen  for  each  of  the  quantities 
n,A,A',A',  and  r.  As  in  Chapter  V  the  numbers  k  are  determined . 
The  expression  for  T   f f  may  then  be  evaluated.  If  the  solutions 
derived  earlier  satisfy  the  energy  equation,  then 


r    =  r 
eff   . 


The  addition  of  a  short  subroutine  to  the  existing 
programs  allows  one  to  perform  these  steps  (see  Appendix  D) . 
The  range  of  r  was  set  from  0.002  to  5.00  and  representative 
values  of  the  other  quantities  used,  over  this  range,  it 
is  found  that 


r 

_§M  »  i 


Thus  the  approximated  form  of  the  energy  equation  and  the 
results  of  Chapter  IV  and  V  are  inconsistent. 


A  Reconsideration  of  the  Approximate  Energy 
Equation 

Non-Static  Solutions  and  the  Adequacy  of  the  Equations 

In  view  of  the  inconsistency  between  the  solutions 
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developed  in  Chapter  IV  and  the  prediction  from  the 
approximate  form  of  the  equation  for  energy  transport, 
the  criterion  for  establishing  this  approximate  energy 
equation  must  be  reexamined.   The  crucial  element  in  the 
choice  of  an  approximate  equation  for  energy  transport  was 
the  identical  satisfaction  of  the  relation  for  the  subscript- 
zero  terms  with  the  substitution  the  static  isothermal 
solutions.  The  analyses  of  the  momentum  and  continuity 
equations  also  resulted  in  equations  for  the  subscript-zero 
terms.   These  equations  were  satisfied  not  only  by  the  static, 
steady-state  axially-symmetric  solutions  but  also  by  the 
non-static,  steady-state,  axial lyrsymme trie  solutions.   Do 
these  non-static  solutions  satisfy  the  subscript- zero 
relation  from  the  energy  equation? 

For  the  non-static  case  the  equation  for  the  subscript- 
zero  terms  is 


v  5  up      -i 

-^PoVuo  ("*  +B-fw)+V^oJ 

Ypovo  .  Bpo  +  fo  .  _!£_  ]      !lzo  =  Q 

"  [- +  —  +  ^T    FNr  J    3z        ' 

2 

where  the  velocity  does  not  contain  a  suppressed  exponential, 
Limber's  non-static  solutions  (see  Chapter  II)  do  not  exist 
in  closed  form,  but  for  any  of  his  solutions, 


m   2tt(1  -  p) 
vuopo  "  u(u  -  p) 
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where  p  is  the  parameter  discussed  in  Chapter  II.  Further 
the  circular  velocity  is  known. 

The  strategy  here  will  be  to  examine  the  next  most 
simple  case  after  the  static  isothermal  envelope — the  non- 
static  isothermal  envelope.  Only  the  region  near  the 
equatorial  plane  will  be  considered.   The  above  equation 
under  these  conditions  leads  immediately  to  the  result 
that  the  quantity 


(422LL  +  B-Jt-  +  *3ij  j2^1  -P>) 

2w         FNr    2      (a)  -  p) 


should  be  constant  within  a  small  second-order  term, 


(  JU£  ?[(2tt)  (1  -  p)  1  m 


VA„      a)  -  p 

(j)0 


However,  the  quantity  above  is  clearly  not  a  constant. 
Therefore,  the  equation  is  not  satisfied  in  this  simple 
case.  The  approximate  form  of  the  energy  equation  must 
be  held  in  doubt. 


Final  Comment 

In  this  chapter  the  distinctive  nature  of  the 
energy  equation  was  discussed.  A  simple  criterion  for 
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choosing  the  approximate  form  of  the  energy  equation  which 
represents  the  same  formal  mathematical  order  of  approximation, 
as  the  analysis  of  the  previous  two  chapters  was  proposed. 
Such  an  equation  was  found.   In  the  end,  however/   this 
equation  was  inconsistent  both  with  the  analysis  of  Chapters 
IV  and  V  and  with  the  non-static,  steady-state,  axially- 
symmetric  solutions. 

The  solutions  developed  in  both  the  previous  steady- 
state  studies  and  in  this  work  depend  on  equations  (continuity 
and  motion)  which  contain  no  radiation  dependent  terms  (save 
for  a  slightly  reduced  stellar  mass  in  the  gravitational  term) . 
This  neglect  has  been  adequately  justified  (Limber  and  Marl- 
borough, 1968) .   That  the  solutions based  on  these  equations 
do  not  satisfy  an  energy  transport  equation  without  radiation 
terms  is  not  an  indication  that  the  solutions  are  invalid, 
but  rather,  that  radiation  must  play  an  important  role  in  the 
energy  equation.   This  is  so  even  though  radiation  effects 
are  of  no  consequence  in  the  other  two  equations.  An  energy 
equation  which  can  be  used  to  improve  the  existing  solutions 
awaits  further  study  of  the  radiation  field  in  the  disk. 


CHAPTER  VII 
•FINAL  COMMENTS 

Stability  of  the  Steady-State  Flows 

The  existing  steady-state  solutions  to  the  hydro- 
dynamical  flow  in  Be  stars  are  approximate  solutions  to 
approximate  forms  of  the  equations  of  continuity  and  motion, 
yet  in  principle  there  is  an  exact  steady-state  solution  to 
the  exact  hydrodynamical  equations  for  any  problem.  With 
regard  to  stars  this  may  be  of  little  value,  for  one  need 
look  no  further  than  our  Sun  to  see  that  a  stellar  surface 
is  not  a  true  steady-state  boundary.   Yet  if  it  were,  and  if 
the  solutions  were  known,  they  would  still  not  necessarily 
exist  in  nature.  There  is  a  more  basic  consideration. 

A  solution  to  the  hydrodynamical  equations  must  be 
stable  against  infinitesimal  perturbations;  a  small  perturba- 
tion, if  it  arises,  must  decrease  in  time.  The  discussion 
of  the  stability  of  the  steady-state  solutions  will  reintro- 
duce most  of  Chapter  IV  in  a  new  light. 

The  complete  procedure  for  the  investigation  of  the 
stability  of  a  fluid  flow  is  outlined  by  Lin  (1967) . 
Generally  only  the  first  stage  of  this  procedure  is  followed, 
and  this  stage  is  formally  the  same  as  the  linearization 
analysis  of  Chapter  IV.   In  short,  the  work  of  Chapter  IV 
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and  V  may  be  used  in  an  investigation  of  the  stability  of 
steady-state  flows.  The  subscript-one  term  may  be  viewed 
now  as  a  perturbation. 

The  general  stability  problem  for  steady  flows  about 
a  finite  body  has  never  been  solved.  On  the  basis  of  the 
special  cases  which  have  been  successfully  treated  a  pattern 
emerges.  The  solutions  to  the  perturbation  equations  are 
those  with  a  time  dependence  in  the  form  exp  (-iut) ,  where 
u  is  a  frequency  which  can  be  determined  only  by  solving  the 
perturbation  equations  with  respect  to  the  prescribed 
conditions.  This  much  is  similar  to  the  analysis  of  the 
approximated  version  of  the  linearized  equations  discussed 
in  Chapter  IV. 

However,  o  being  a  real  number  is  generally  the 
exception  rather  than  the  rule.   Depending  upon  the  problem, 
o  may  have  imaginary  parts  which  lead  to  damped  perturba- 
tions or  positive  imaginary  parts  leading  to  solutions  which 
increase  exponentially  with  time.  These  last,  of  course, 
are  the  hallmark  of  an  unstable  flow.   Insofar  as  the  steps 
leading  from  linearized  equations  to  the  k  equations  are 
valid,  the  steady-state  flows  are  neither  absolutely  stable 
nor  absolutely  unstable. 

The  flow,  though  the  steady  state  component  is  not 
absolutely  unstable,  can  no  longer  be  regarded  as  a  steady- 
state  flow.  While  the  solutions  of  Chapter  IV  were  not 
fully  determined,  they  were  not  required  to  have  infinitesimal 
amplitudes — nor  is  there  any  reason  to  suppose  that  they  will, 
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The  general  solution  for  each  flow  variable  is  still  the 
sum  of  the  steady-state  term  and  the  perturbation.  One 
has  passed  into  the  region  of  periodic  flows-  (in  the 
sense  that  the  subscript-one  flows  are  periodic) ,and  the 
steady-state  flows  are  at  best  approximations  (though  they 
may  be  good  ones) . 

Logically/  one  must  then  proceed  to  the  question  of 
stability  of  even  periodic  flows.   In  principle  one  could 
perform  a  stability  analysis  on  the  subscript-one  terms  in 
much  the  same  manner  that  was  followed  in  the  analysis  of 
the  stability  of  steady-state  flows.   Such  an  analysis  is 
characterized  by  exceptional  mathematical  difficulty,  and  is 
seldom  attempted  for  even  the  simplest  flows.  None  will  be 
attempted  here.   If  the  periodic  flow  were  unstable,  the 
flow  would  be  turbulent.  A  treatment  of  turbulence  is 
beyond  the  scope  of  this  work. 

Conclusions 

It  was  concluded  that  the  motion  and  distribution 
of  the  circumstellar  matter  about  typical  Be  stars  could  be 
treated  hydrodynamically .  Typical  disk  conditions  placed 
the  flow  within  the  realm  of  hydrodynamics,  but  close  to 
the  limits  of  validity  of  that  discipline.   It  was  also 
shown  that  over  most  of  the  disk  Euler's  equation  could  be 
used  for  the  equation  of  motion. 

Each  flow  variable  was  written  as  the  sum  of  a  known 
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steady-state  axially-syrometric  term  and  an  unknown  term 
containing  the  temporal  and  angular  dependence.   The  squares 
of  these  latter  quantities  were  assumed  small.   Using  this, 
the  equations  of  motion  and  continuity  were  reduced  to  a 
set  of  linear  equations  for  the  terms  with  temporal  and 
angular  dependence.  A  static  isothermal  solution  from  the 
literature  was  chosen  for  the  known  steady-state  axially- 
symmetric  term.   Solutions  were  found  to  a  somewhat 
approximated  version  of  the  linearized  equations.   In  this 
process  two  parameters  were  introduced.  While  the  solutions 
could  not  be  fully  determined,  the  dependence  on  both  time 
and  angle  is  explicit. 

The  angular  dependence  entered  the  solutions  in  an 
extremely  simple  form  involving  an  integer  n,  but  understanding 
the  temporal  dependence  required  computation  of  the  quantity 
called   | PER |  ,  which  represents  the  periodicity  of  the 
time  variation.  These  machine  computations  were  carried  out 
for  a  wide  range  of  parameter  values.   The  results  were  that 
these  quantities  (|per|)  were,  under  most  conditions,  real. 
Further,  their  dependence  on  n  and  oj  was  generally 


3/2 

a)  : 


These  calculations  were  in  qualitative  agreement  with 
relevant  observations. 

The  third  equation  of  hydrodynamics,  the  energy 
equation,  was  examined  insofar  as  could  be  in  the  absence 
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of  an  effective  treatment  of  the  radiation-fluid  interaction. 
This  approach  was  found  to  be  inconsistent  with  earlier 
work.   The  conclusion  of  the  matter  was  that,  in  the  absence 
of  some  accommodation  of  radiation  effects,  no  equation  for 
energy  transport  was  available  which  could  improve  the 
existing  solutions. 

Finally,  the  steady-state  flows  were  examined  for 
stability  effects.  While  the  steady-state  flows  were  not 
found  to  be  absolutely  unstable,  the  most  general  solution 
appeared  to  be  periodic.   The  stability  of  periodic  flows 
was  not  examined. 


APPENDICES 


APPENDIX  A 
Exchange  of  a  Differential  and  an  integral  Operator 


Let  f,  (w)  and  f,^)  be  functions  of  w  such  that 


f^w)  =  -  f2(w)  ; 


further,  let  g(u,<S>,z)   be  a  well  defined  function  in  the 
interval  from 


z  =  f^tw) 


to 


z  =  f2^ 


What  is  to  be  shown  is  that 


f2(w)   a  a    f2<w) 

L    ,  x  3^(u'*'z)  dz  =  ^  f     g(u>,4>,z)  dz 


Define  the  quantity 

z' 

H(u),<J>,z')  =  /  g(w,<|>,z)  dz  +  constant. 


Now 

f2(w)  f2(w)  ^(w) 

/       g(w,<J>,z)  dz  ■  /      g(w,<j>,z)  dz  -  J   ~   g(w,<j>,z)  dz 
f  !<<*>) 
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=  .H(W,<j>,f2)  -  HM,*.,^)   ; 


hence 


3   r  f2^  •  a 

^  J  g(w,<f>,z)  dz  =  |^  [H(w,<j>,f2)  -  H((a,4>,fif)] 


In  the  partial  differentiation  the  to  dependence  that  enters 
through  the  functions  f^w)  and  f 2  (cu)  can  be  formally 
separated  from  the  general  u  dependence;  this  is  shown  below 


bw.M.f,)--(-^f1)I-^^-1^^ 


+  ff*«,*,f2>  -  ||H(a,f*,  r )    f 


where  the  bar  over  f±   and  f2  symbolizes  that  the  partial 
differentiation  treats  these  as  unrelated  variables. 
But 


3        —     3H      —     a     9 

*1 


and,  since  the  partial  differentiation  does  not  "see"  f 


and  f2  , 


j\    £^  f  _ 

5tt  /   g<wr*,z>  dz  =  /   ^|  g(u>,«J,,z)  dz 

1  *1 
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It  remains,  then,  only  to  show  that 


9f      '  3f 


Since 


f-^co)  =  -f2(w)  / 


it  follows  that 


and 


Clearly 


hence 


and 


3f  3f 

^r  c-i  =  -  to2-  (w) 


9  _3_  3f  2  _3_ 

9^      9f2  3^  "  "  3f, 


9H    9fl  3H      3f2 


9H  ^£2        _9H      *£l        n 

TIJ     9oj  "    3^        3u        U 


w~r  /      "  g(w,<|>,z)    dz  =  /  4_:g(w,<J>,z)    dz 

3W  f^U))  f^U?1 


APPENDIX  B 
Separation  of  the  Equations  and  the  ic  Relation 

The  three  equations  which  contain  the  dependent 
variables  are 


Luwl  "  2fV  -  ^7  ai  ' 


f         C2   3al 

Lu4>l  +  2  ua)l  =  —  Jf 


La,      +  -^+  i.**i  =   o  . 

1     w    w  9<J>         ' 


In  order  to  apply  traditional  mathematical  procedures  to 
this  problem,  it  is  useful  to  construct  separate  equations 
for  each  flow  variable. 

Using  the  third  equation,  one  can  write  that 

%i  ■  -4Lcri  +  £af  u*i]  •' 

This  may  be  used  in  the  remaining  two  equations  giving 
[~L  4"2f]    «♦!      =    I>2  +  fJ    ]«i 


151 


152 


and 

.  c. 


r .  _  .  •  f  '  3  "I       r  fwL  •'  _2   3_  1 

L  L"2.3*J^rL"2"+u  3<p  J    ai 


The  operators  which  precede  a,  in  each  equation  are 
allowed  to  operate  on  the  other  equation.   Since  these  two 
operators  commute,  the  right-hand  sides  of  the  resulting 
equations  are  equal.   Subtracting  one  from  the  other  yields 

c2L  32     2c2f   B    g^fL2.  _3  _   2   ,   _ 
r  T3    fwL2   3   ■  Cl  T   •  Clf   3  S  i  „    _   n 


or 


,    c0L   „2       c      ,        2c9f   c  f   . 
(  .L3   ♦  -V-  V,f  (■  -4  +  *f2  )  L  +  (-4-  -  -Ij)^ -,.}u  -  0 

3  cp     2w  2w 


Let 


u<bl  =  Unk  exp  ^in*  +  ik(w»n)t]  * 

and  substitute  this  into  the  equation  above.  The  result 
of  these  steps  is 

i  {  -a)(k~  +  nf)3  -  -2—  (k  +  nf )  +  (  -r  +  *»f2  )  (k  +  nf) 

c. f   2c0f 
_  (1  +  — ^-  In}-  0  . 
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Hence, 


tok3  +  3confk  +    (  —  -   tof2  +  3wn2f2  -  -|  )    k 

CO  CO 


3 
c_n  f        2c.„fn  ,  _   _     .  c.nf 

+    (  -?■ 1 tof3n  +  ton3f3  -  -±r-  )      =  0 

1   w       w  2co4 


Much  of  this  analysis  may  be  used  for  the  reduction 
to  an  equation  for  g,.  The  operators  which  precede  u  ,.  in 
each  of  the  two  equations 

(  -LfdT  "2f  )  »*i<!  (wl2  +  ~l)  al 


and 


I  L"  2  8?  J    U(f»l  -  I  T"  +  ~  3?  J  CT1 

are  allowed  to  operate  on  the  other  equation.   Since  the 
operators  commute,  the  left-hand  sides  are  now  equal.   Subtrac- 
tion and  some  rearrangement  gives 

.j        c0L     ~2  c.  0  2c0f       c,f,   , 

-  ^3  + .4-  rrf +  ( r* +  Mf  J L  +t-*- "  ;M  I? }  °i "  •• 

dtp  2d)  ZU) 


The  steps  to  the  k  equation  need  not  be  repeated. 

The  third  of  the  three  original  equations  can  be 
rewritten 

Lai ~  "  »  57  u<pi  • 
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Substitution  into  the  remaining  equations  gives 


2   Cl  V„   _  r   ^T   °1  9 


(  L  +  -T  Kl  "  (  2fL  "  -T  k  )  V 

0)  0)     r 


and 


0)   0<p 


The  mutually-commuting  operators  preceding  Ux,  in  each  equation 
are  allowed  to  operate  on  the  other  euqation.  On  subtraction 
and  rearrangement  one  has 

,        c_L      a2  c,  0  2c_f        c.  f        «  u  , 

{wL3  +  _2_     3      +    f        1  +  wf2   j    L  +r  _2_         1         ij    (jol)    „  0 

w      9<j>2  2co4  w  2co4      ^  u 


APPENDIX  C 

Evaluation,  of  the  Determinant  Formed  from  the  Equations 
Relatxng  5^7  C^,  S.^ 


The  determinant  which  must  be  evaluated  is 


i(k  +  nf) 


-2f 


-c^/ftf 


f/2 


i(k  +  nf) 


-inc2/w 


1/0) 


in/w 


i(k  +  nf) 


On  expansion  in  terms  of  2  x  2  determinant,  one  has 

i(k  +  nf) 


i(k  +  nf) 


in/w 


-inc2/w 


i(k  +  nf) 


+2f 

f/2 

1/U) 

-inc2/w 
i(k  +  nf) 

cl 

"7 

This  gi 

Lves 

f/2 


1/0) 


i(k  +  nf) 


in/w 


n  c->         -  f  ,<r   _^v  .  n°2 


i(k  +  nf )  (  -(k  +  nf)2 f-  )    +  2if  (  §  (k  +  nf )  +  -^  ) 
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4  .  '    V  2(ii  to  J 

CO. 


or 


2 

i    {   k"3+   3nfJ?2  +    (  SS-  -  f2  +   3f2n2   -  ^l  )      Jc 

to  u 


c  n3f        2c   fn  3  -   3       c^f 

+  ( .-^-2 2~  "  f  n  +  f  n     -  — y  J    1    • 


The  quantity  in  brackets  is  the  equation  for  k.  The  determinant 
is  identically  zero. 

If  atfion-zero  2x2  determinant  can  be  formed  from 
matrix  elements  above,  the  matric  is  of  rank  2.   There  are 
three  2x2  determinants  which  appear  in  the  expansion  of 
the  3x3  determinant  above.   The  condition  that  all  three  be 
zero  is  that 

2.2 

c2  =  -  oj  f 

Over  most  of  the  disk  such  a  value  represents  a  largervalue 
of  T(-c9)   than  those  in  use  here.   For  the  computations 
described  in  Chapter  V  the  matrix  is  of.  rank  2. 


APPENDIX  D 
The  Program  CALSOL 

The  major  variable  names  for  the  program  are  related 
to  the  symbols  used  in  this  work  in  Table  3.   The  program 
has  four  major  divisions: 

(1)  Main  program, 

(2)  Subroutine  POLRT, 

(3)  Subroutine  ALPl, 

(4)  Subroutine  BETl. 

The  second  division  POLRT  is  essentially  unchanged  from  the 
discussion  given  in  the  IBM  publication  System/360  Scientific 
Subroutine  Package  Version  III  -  Programmer's  Manual.  The 
main  program— minus  COMMON  and  DIMENSION  statements— is  given 

in  Table  4. 

There  are  internal  control  statements  in  the  program 
which  change  from  run  to  run.  When  these  occur  in  the  program 
the  expression  control  'x'  appears  where  'x'  is  a  number 
The  locations  of  these  statements  is  marked  by  a  double 
asterisk  (**)  to  the  left  of  the  statement.   CONTROL  1  and 
CONTROL  2  together  control  the  values  or  sets  of  values 
which  c,  and  c,  are  to  assume.   CONTROL  3  tells  the  program 
at  what  values  of  n  the  calculation  is  to  be  performed,  and 
CONTROL  4  chooses  the  points  in  the  disk  for  which  the 

computations  are  to  be  performed. 
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Subroutine  ALP1. 

This  program  calculates  the  quantity  r  --  discussed 
in  Chapter  VII.   The  program  is  shown  in  Table  5. 

Subroutine  BETI 

This  subroutine  calculates  the  density,  the  circular 
velocity,  and  the  radial  velocity  in  that  order  when  given 
S  ;-.   This  program  is  shown  in  Table  6.   Again  all  secondary 
statements  have  been  deleted. 
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TABLE  3 
The  Program  Names  for  Important  Quantites 


SMAS 

M, 

SRAD 

= 

R 

< 

SVEL 

= 

V 

i 

RMU 

= 

V 

RTEMP 

= 

1 

< 

PvLAM 

= 

A 

RGAM 

= 

r 

Q 

= 

Q 

PSI 

= 

* 

CI 

= 

c 

s 


s 


C2 

)B 

C2 

OMEG 

= 

0) 

N 

= 

n 

ROYTR 

= 

Real 

part  of  k 
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TABLE  6 
Subroutine  BET1 


ARG1=(6,283*10.00)/PER 

ARG2= (3. 1416/4. 000  )*I 

SK (J.N. )=1. 000/(1. 000+j/ (10. 00)) 

DO  80  IPHC=1,8 

DO  90  ITIC=1,10 

ITIM=ITL)-1 

IPHI=IPHO—l 

DENN (N, IPHO , ITIO) =SK ( J,N) * (COS ( (ARGl*ITIM) + (ARG2*IPHI) ) ) 

SOLN(N,IPH0fITI0)=(SK(J,N)*Y)*(COS(ARGl*ITIM+ARG2*IPHI)) 

RADN (N, IPHO , ITIO) =- (SK ( J, N) *XBAR) * (SIN (ARG1*ITIM+ARG2*IPHI) ) 

90  CONTINUE 

80  CONTINUE 

RETURN 

END 


APPENDIX  E 
Pertinent  Integrals 
1.   The  integral 


-z_  20) 


4   4 
where  z  /o)  and  higher  powers  of  z/co  have  been  neglected. 

For  the  region  of  interest  (|z|/o)  <_  0.1)  these  may  be 

neglected  even  under  the  integral.   Using  the  definitions 

of  Chapter  IV, 

+zo 

r  1  3pl  -   .  f   1    A   ,   3al 

/    ?3Tdz  -  I  5  "  73  J   3t 

-z  2o) 

o 


2.   The  integral 


+z  +z    a 

LI     "^r  fev^ijdz  -  -    /     (|  .-^p  +  ?  vwl)  dz 

—  z  —z 

o  o 


i    /    a.  v  .  dz  =  i   /    (1  -  S-*)^"  <az  + 
o  o 
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i      "^  o  +z        2 

-    r   °  l  r,  _  z_)  v  .  dz  -  i   f   °:4lA)  vMldz. 

0)    J      JJ  l1      2^  vwl      to    /     «<•   2^2^   Oil 


Again,  the  terms  in  z  /w  are  so  small  that  they  can  be 
safely  neglected.   This  expression  reduces  to 

+zo 

ST     /    3£   S  V0)l)  dz  "  l» :J3S-  +  H4  Vl' 

2co 
-z 
o 


3.   The  integral 


+z 
+z  o 

/    rw  h  V4l   dZ  =  I?  _z  /   riV  dz 
~zo  o 


+z        ,      2 
3       o  v.,  r   1    z  i  - 

to  terms  in  z  /w  and  higher.   This  expression  gives 


J    ra3   3<J>  az  lu>        2aJ4^   5? 


4.  The  integral 

)  o 

to  the  usual  approximation.   This  can  be  written 


+z„  +z«  9 

F3*  Pldz-3*.z/     Plfe"^ 
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-■  J     r  I?  dZ    U   ^  3* 


5.   The  integral 

+  2 


-z  ;     3z  *  a)  ' 


reduces  to 


2vzl<zo> 
(I) 


immediately. 


APPENDIX  F 


An  Equation  for  f  ff 


The  quantities 

uwl  =   2     Unk  exp    [in<|,+  i5?t] ' 
n,k 

U<J>1  =  E     Cnk  exp    tint+ikt]  , 
n,k 

and  ff1     =  I     Sn^  exp    [in<J)+ikt] 

n,k 

from  Chapter  IV  are  substituted  into 

uf  3t  V  +^-T-  +  B      ¥m  +  =T3  +  2  J   at 


3u, 


■frt£^?     Wl"   "^  5        *        S?        2FN.4J    3* 


(-T  +  Bf        FNS  +   2w3FN  2   J    3* 


An  operation  of  the  type 

+T  +T" 


Lim  1 


—       _/     dt        /     d(j>  exp[sin<|>-ikt] 


T   -»»  oo      4ttT      -T  -it 
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is  performed  on  the  equation,  and  integrals  and  sums  are 
interchanged.   Finally  the  equation  is  divided  by  Snk, 
and  the  ratios 


.-   Unk 
bnk 


Cnk 


Y  "  SnZx 


introduced.   The  result  of  these  steps  is 


*  ♦  (t1^  +  B  +  d^r  -  rfc +  *&  k  - 

2FNd) 

rB       5^  2T'  i  -       r3tof2       B       5<J>  _  _l .  _A!_)     nv 

-    rw2f3        „*  f  Af       _5f£i    n 


The  quantity  r  where  it  appears  explicitly  in  this  equation 
is  called  r&^f .   Solving  for  r#ff  one  finds 

2    2  o. 


r_  «  {-  (  M"  +  *  fc  o!  +  TT    * 

0)        Qw 


3oj3  '  n       i«      25      w4  J 


1      2  (0        2C0        Qw2        JT  ^      y 
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2  u  3w 


r 3k        5fn  •) 
A   2  +      2     ' 


where 


and 


'3         Q   A 


=4=!  a-. 
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